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THE ZETA FUNCTIONS ARISING FROM 
QUADRATIC FORMS, AND THEIR 
FUNCTIONAL EQUATIONS 
By L. J. MORDELL 
[Received 7th March 1930] 


§ 1. The functions of a complex variable s defined by a Dirichlet’s 
series in which the general term involves ¢-*, where ¢ is a homo- 
geneous quadratic form in the variables of summation having for its 
real part a definite form, are well known in analysis.* For one 
illustration, take 
d(x, y) = A(x—w, y)(w@—wy y) = Ax*+ Bay+ Cy’, (1) 
where w,, w, are any two unequal complex numbers neither of which 
is real, and the notation is such that J(w.) < I(w,), or, if these are 
equal, R(w,) << R(w,). The discriminant of 4(x, y) is given by 
A= 4AC— B? = —A?(w,—w,)*. (2) 
The multiple-valued functions are defined by 


£8 — es log ifi+is ares —am <argé <7, 


but, if € is a real negative number, its argument +7 will be stated. 


Again, (a, y)® = AS(x—w, y)*(x—wy y)’, 
so that arg 4(x, y) is defined as the sum of the arguments of its three 
factors. 

In particular, if 2 and y are real, then, when y > 0, 

(x—w, y)? = Y*(x/y—a)*, 
where 0 > arg(x/y—a@,) > —7, if I[(w,) > 0, 
and 0 < arg(2/y—a,) <7, if I(w,) <9. 
When y <0 pn | elie . 
"7 @—eyy)? = lyl*(@/ly|+-e), 
so that now 
0 < arg(x/y+w,) <2, if I(w,) > 0, ete. 

When y = 0, d(x, 0) = A*(x?)s, 
the argument of 2? being zero. 

Two essentially different cases arise according to the signs of 

* Some references are given in my paper ‘Kronecker’s fundamental limit- 
formula in the theory of numbers and elliptic functions and similar theorems’, 
Proc. Royal Soc., 125 (1929), 262-76, which is referred to as R. 8. 
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78 L. J. MORDELL 
I(w,), I(we). These signs may be different or the same. There is no 
loss of generality in supposing that I(w,) >0. Denote by vA that 
value of the radical such that R(A/vA) and J(w,) have opposite signs. 
Suppose first that I(w,) >0, J(w.) <0, which is apparently the 
only case considered by other writers. Put 
VA (x, y) = Ax?+ Bay+ Cy’, 
then the real part of (x,y) is a positive definite form. For, since 
here R(A/vA) > 0, (2) gives, on taking y = 1, 
(a, 1) = is aye (3) 
Write w, = a,+7b,, wo = a,+%b, and so 
—i[a,—a,+1(b, —b,) (x, 1) 
= (w—a,)(~w—a,) —b,b, —i[b,(~—a,)+b,(%—a,)]. 
On ignoring a positive denominator, the real part of (2,1) is 
(b, —b,)(x—a,)(x—ag)+ (a, —a)[b,(~—ay)+b,(%—a,)]—b,b,(b, —b,). 
If we write = x—a,, n = X—Qg, this is 
(b, —by)En+(n—E€)(byn +b, €)—b,b.(b, —b2) = bn? —b, €&?—bbo(b, —,), 
and the result follows, since b, > 0, = <0; 
This means that 


—tn <argy(x,y) < 


and so there is no need to factorize We, ) in order to find arg (x, y) 
as already defined. For 


—thn <arg —i(w,;—w,) < $2, 


and a figure shows that 

—n7 <arg(w—w, y)+arg(x—wyy) <7. 
Hence arg x(x, y): ), defined as the sum of the arguments of its factors, 
lies between —3z and 37, and since ja ” has a positive real part, 


—tr <argy(x < 4x. 


f(s, A, B, C)=f(s) . >” d(m,n)-, (4) 


n oo m= nH 


Put 


where m = n= 0 is excluded from the summation. 

The series converges absolutely for R(s) > 1, and the function of 
s so defined is really analytic throughout the finite s-plane except 
for a simple pole at s=1. It satisfies the functional equation 
I'(1—s) 


fl) = (2x VA 


f(l—s), 
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where R(A/VA)>0. The proof yan here shows incidentally that 
for the usual zeta function {(s) = > n= 

8) = 7®-t — = __ = {(1—2). 5.1 
(s) fay (5.1) 


It is shown that is tantamount to the relation 


(5) 
rf 2Qm7rix dx 
(Aa?+ TBs n+ Cn?)s 


e2naix dx 
6 
| (Aa?+ Bam+ Cm?)!-8’ (6) 


where m, » are both real and different from zero, 0 < R(s) <1 and 
(Aa?+ Bay+ Cy*)/vA 
has for its real part a positive definite form. The last condition is 
essential in applying the idea, given later, which obviously suggests 
itself when we try to prove (6). 
Take next the case [(w,) > 0, [(w.) > 0. 


oo 


Put fils) =f,(s,4,B,C)= S > dlm,n)-*, (7) 


m=-0 n=1 
for it will be seen that these limits of summation in are more 
convenient than +00. 
The series converges obeclatally if R(s)>1. For put 
d(m,n) = aym?+-bymn+c,n?+7(aym*+-bymn-+-cn*) 
= $,(m,n)+ip,(m, n), 

where @,,5,,@,... are all real and ¢,(m, ”), (m,n) are both indefinite 
forms. Then, for all integer values of m, n other than m= n= 0, 
d(m,n)| has a positive lower bound* and (m*-+-n*)/|(m, n)| 
is bounded. For let 6, be a root of ¢,(x,1)—0 and @, a root of 
d(x, 1)=0. There is no need to consider the terms arising from 
n = 0, since > m-** converges. Then, given « > 0, consider separately 
the terms arising from |m/n—8@,| << «€ and |m/n—6,| >. For the first 
set, d,(m, n) may be small, but since the roots of ¢(x, 1) are complex, 
d(x, 1) = 0 and ¢,(x, 1) = 0 have no common root. Hence ¢,(m, n) 
will not be small, if the notation is such that 6, is the greatest of the 
four roots of ¢,(x, 1) = 0, ¢.(z, 1) = 0 and « is sufficiently small; and 
then clearly lbo(m, n)| > K(m?+-n?) 


* By positive lower bound we shall understand essentially a bound greater 


than zero. 


G2 





80 L. J. MORDELL 
for a suitable constant K. For the terms given by |m/n—0,| > we 
obviously have 
lp, (m, n)| > K,(m?+n?), ldo(m, n)| > K,(m?+-n?) 
with suitable constants K,, K,. Hence the result. 

But now an entirely new function arises. The real part of 
d(m, n)/VA is an indefinite form. The definition of the arguments is 
such that arg 4(m, n) does not lie between -+-7, since, for given y > 0, 
arg(v—w,y) and arg(wz—w,y) both decrease from 0 to —z as x 
decreases from 00 to —oo. I have shown that /,(s) is really an integral 
function of s, for it will be seen that in R. 8. terms have been included 
arising from n = 0 which do not arise here. 

§ 2. The question at once suggests itself: Does f,(s) possess a func- 
tional equation of the type (5)? The usual methods cannot be 
successful for they depend upon theta functions, wherein it is essential 
that the product of ¢(m, ) and some constant should have for its 


real part a positive definite form in m,n. Recently* I have developed 


another method wherein the theta function is replaced by an applica- 
tion of Poisson’s summation-formula. Then on utilizing an idea due 
to Siegelt I have found{ with comparative ease the functional 
equations of all the zeta functions in an algebraic field. In principle, 
the method is applicable when the old methods must fail. 

But, applied so far, it depends upon a preliminary result of the 
type (6), wherein, for the usual proof, the definiteness of the quad- 
ratic form is essential. I have, however, found another proof of (6) 
which also gives the corresponding result when I(w,) > 0, [(w.) > 0, 
and hence also the functional equation for f,(s). It is not so simple 
as before and now takes the form, if 0< R(s)<1, I(w,) > I(ws), 
R(w.) > R(w,) when I(w,) = I(w,), 


; o..1 1 (1—8), 
fils) = 2nvay Te" f,—s)— 


—. | | e2zmnrix dx 


— (7.1 
(Ax2+Bre+oy ‘ ) 


— 27 sin 78 cos 378 —s 
a 
n=1lmel 


‘ 


w1 

* *Poisson’s summation formula and the Riemann zeta function’, Journal 
London Math. Soc., 4 (1929), 285-91. 

+ ‘Neuer Beweis fiir die Funktionalgleichung der Dedekindschen Zetafunk- 
tion’, Géttinger Nachrichten, 1922. 

t ‘On Hecke’s modular functions, zeta functions, and some other analytic 
functions in the theory of numbers’, Proc. London Math. Soc., (Read Dec. 12, 
1929; to appear shortly). 
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where the path of integration is a straight line parallel to the positive 
direction of the real axis of x, or slightly inclined to it at a positive 
angle if I(w,) = I(w,). For other values of s, we need only replace 
the linear path of integration w,, 00 by the usual cut going from oo 
round w, in the clockwise direction and returning to 00, and divide 
the new integral by e?7*—1. On the lower side of the cut at x = 0, 
when J(w.) + I(w,), we can take 


arg(Ax?+ Br+C)=arg A. 


We first prove (5). As (x, y) can be used instead of d(x, y), we 
may suppose that A= 1 and that ¢(x, y) has for its real part a 
positive definite form. In (4) write separately the terms given by 
n = 0 which are 

2 > A-“m-* = 2A-*f(2s). 
m=1 
For given n+ 0 apply Poisson’s summation-formula to the summa- 
tion in m, whence as in R. 8. 


ae - e2mriz dy 
f(s) = 2A~*C(28) + > >. | Azx?+ Ban+Cn?)s’ 


=-o m=-0 ©, 


the path of integration being the real axis of x. 
Isolate the terms with m=O and write nx for x in the corre- 
sponding integrals. Then 


a) 


py I(m,n,8), (8) 


f(s) =2A-*(2s)+2¢(2s—1) )| + 
ae BO y m,n=—@ 


4, 


where both m = 0 and n = 0 are excluded from the summation and 


I(m, n, 8) = (9) 


e2mrizx dy 
(Aa?+ Ban+ Cn?" 


By writing nx for x and integrating repeatedly by parts, it is easy 
to show that J(m, n, s) and >’ I(m, n, s), defined for R(s) > 1, are 
really integral functions of s. Since the behaviour of f(s) is known, 
(8) enables us to study f(s), as is done in R. 8. 

There is no need, however, to assume any properties of {(s), since 
they are implicitly contained in (8). Take B=0, A>0, C>0, 
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interchange A and C, and subtract. Since f(s) is symmetrical in A 
and C and 


> 


. de nt T(s—}) 
| aero ~ Alcs T(s) 
then i 
(A-*—O-*)£(28) + (4410-4 CA) PY) P(8—3)E(28—1)/T(8) 
= G(s), . (10) 
where G(s) is an integral function of s. This result is really a recur- 
rence-formula between {(s), ¢(s—1) from which it is easily shown 
that {(s), defined for R(s)>1, is analytic throughout the finite 
s-plane except for a simple pole at s=1. Thus, if R(s)>4, (10) 
shows that (As-1— C-1)£(28— 1) 7 (s—4)/T(s) 
is analytic, and hence 
(Al6-D— CHO-D)C(s)P(48)/P(h8-+-4) 
is analytic for R(s) > 0. Since A, C are arbitrary, f(s) is analytic for 
R(s) >0 except possibly for a simple pole at s=1. That s=1 is 
actually a pole follows from Landau’s theorem on the singularities 
of a Dirichlet’s series with positive terms. The functional equation 
(5.1) then shows that ¢(s) is regular throughout the rest of the s-plane. 
Suppose R(s) > 1, then, by de la Vallée Poussin’s theorem* since 
the double integral converges absolutely and can be evaluated by 
integrating in either order, 


» 8 
I(m, n. s) a 27) | | emmn 2rt( 427+ Bnz+e "n”)4s -1 didx. ( 11) 


( 
rs) 


—o 0 


Integrate first for x and note that 


r in a 214 
feve-e de =F ety, R(p) > 0, R(p') > 


e-2ntCn® 2 4+-a2(mi—tBn)* 2Alps—3 dt. R(A?) >0 


fi am? natn® mBmn\,. , 
| exe(— Fa T one )e i dt (11.1) 


* Cours d’ Analyse Infinitésimale (2), 1922-23. 
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As m #0, n #0, the right-hand side is really an integral function 
of s. Write 1/t for ¢ and this gives (6) at once with A= 1, namely 
I(s)I(m, n, 8) = (27)*4T(1—s)I(m, n, 1—s). 
This result, proved for R(s)>1, holds for all values of s, since 


both sides represent integral functions of s. Then (8) gives at once 
that, if 0< R(s) <1, 


i. 
Ax?+- Bu+C)s 


f(s) -24-*(28)— 2424-1) | 


= (2778-1 tor . —2A-1-9f(2— 298)— 


nD 


st aster :|: 
or say = 
f,.—T,—T,=1T,—T,—T,. 
Clearly 7',— 7’, is a symmetric function of A and C. Also (7',—T7',)A’, 
(T,—T;)A!-* are both independent of A and C, since, from (11.1), 
A1-*I(0, n, 8) is independent of A and C for R(s) > 4 and so also, 
by analytic continuation, for 0 << R(s) <1. Hence 
T,—T,=0=T,—T,=T;—T;. 
The first of these is (5) with A= 1 and 0 < R(s) <1. It holds for 
all values of s, since, on dividing by I'(1—s), both sides of (5) are 
integral functions of s. 
The other equations give the well-known functional equation for 
the ¢-function. Thus from 7',— 7, = 0, on putting A= C= 4, B=0 
and 4s for s, we have 


ls) (202 [ de 


1-—is 


and this is (5.1). 

These results proved for 0 < R(s) <1 subsist for all values of s. 

§ 3. Now suppose this method of proof is applied to the func- 
tion f(s). Equation (11) cannot be used, since it is essential that 
Ax?+ Bxry+ Cy? should have for its real part a positive definite form. 
Clearly, when (6) is expressed in terms of w,, ws, the real axis of w, 
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is a line of essential singularities of the integrals considered as 
functions of w,. Hence the result must be found in some other 
way. By deforming the path of integration in (6) into an infinite 
semicircle, and drawing cuts from nw,, nw, to infinity parallel to 
the real x-axis, and by writing x+-nw, for x, etc., integrals arise of 
the form 


bre 
ee da 
x*(ca--b)s’ 
0 
and a similar form with s replaced by 1—s. The relation between 
them is derived from the following theorem given by Dirichlet* when 
c= 1. Suppose a, b, c, p, q have positive real parts, then 
2 


e-buyi-1 dy = P ear yp -1 dx 
aces om Fg) | (13) 
(cya)! (ca-+-b)4 


0 
where 
—4n <arg(cy+a) < 42, —4hn <arg(cv+b) <4n. 


It is deduced immediately, if we evaluate 


co fm 


I e~ery—az—byyp—lya-1 dady 


0 0 
by integrating first for 2 and then first for y. For the double integral 
converges absolutely and de la Vallée Poussin’s theorem justifies the 
interchange in the order of integration. Take p= 1—s, g=s8 where 
0 < R(s) <1, whence 


io 2) 


= e—be dx 
3. I (1—s) E “6(ca-+-a)!-8" (14) 


eax dy 


x(cu-+-b) 


This result must be extended to the case when a and b are purely 
imaginary and R(c) < 0 as well as R(c) > 0. Considerations of con- 
tinuity and analytic continuation suffice for this. Put a= —im+A 
where m +0 and is real, A>0O. Then both integrals in (14) are 
continuous functions of A for A >0 and so we may take A = 0. 
Thus for the right-hand side 


(ca-+-a)s-! = (cx—im--A)§-1 = (cwu—im)§-!+ A.O|(cx—im) |§-? 


* ‘Sur les intégrales Eulériennes’, Journal fiir Math. 15 (1835), 258-63, 
or Werke, 1 (1889), 278. 
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where the constant implied in O depends only on s for small A, since 
cx—im| has a positive lower bound, if R(c) >8 >0. Clearly 


rs —R(b)x 
A | . dx >0asA—->0, 


x1-§|\ca—im|?-* 
0 


since the integral converges if R(b) >0. On the left-hand side, write 


= f+f. 


0 
1 1 
eimr—Ax dx eimr dx 
7 s —- 8 | 3? 
x(ca+-b)' xs(ca+-b)s 
0 0 


since the integral is seen to converge uniformly for R(A)>0 on 


io a) 


0 


comparison with 
1 


dx 
axs(cx-+-b)s" 
0 
The result also follows on writing 
e-4r — 1— Axe-PAx where 0<6< 1, 


since 
1 


x8(ca+-b)§ 
0 0 


ca+b)s| 


1 
imx—O Ary dy 
ali* ir | <4] Bind >0asA>0. 
|ars( 


oo 


eimx—Axr dx al * d(e*m2—42) 
~— | as(ca+-b)*(im—A)’ 


x*(ca+-b)s 
1 1 


Integrate twice by parts. Three integrals arise of which one is 


ie 4) 

eimx—Ar dy 
| axs+2(ca-+-b)s 
1 


which is seen to converge uniformly for A >0 on comparison with 
wo 
dx 
as+2(ca+-b)** 
1 
Alternatively, we can use (14.1). Similarly for the other two 


integrals. 
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This proves (14) with a= —im. The work shows that we may 
now put in this b = —in with n real and not zero. We have then 


¥ elmer dx * eine dx 
I'(s —_ == ]‘(l—8s : , 15 
(8) , ( | aera es 
0 0 


if m, n are real and not zero, R(c) > 0, and 


—4n < arg(cx—in) < 3a, 
where, when x = 0, 


arg(—in) = lim arg(cx—in) = 


z—>o0 
according as n <<0orn>090. Also 
—4n < arg(cx—im) < }r. 


The first integral in (15) defines an analytic function of c, not only 
for R(c) > 0 but for all points in the c-plane slit along the imaginary 
axis from 0 to +700 according as n > 0 or n <0. Considerations of 
continuity show as before that, if R(c) = 0 and nJ(c) < 0, then 

arg(cyv—in) = +47 
according as n <0 or n>0. The theory of analytic continuation 
then shows that, when we remove the restriction R(c) > 0 and 
exclude the part of the imaginary c-axis given by nJ(c) > 0, we 
must have 
hor, if n> 0, 


° 


—3n <arg(cx—in) < 


—4n < arg(cx—in) < 3n, ifn <0. 


The same remarks apply to the second integral in (15) and hence 
if m, n are real and not zero, 0 < R(s) <1, then 


D - 


: eime dx . { eine dx . 
I‘(s) : —— = [(1—s) — : (16) 
x8(cxu—tn)® x1-*(cx—im)}-8 
0 0 
for all values of c excluding the part or parts of the imaginary c-axis 
from 0 to imoo and 0 to inoo, provided the argument of cv—in lies 
between —3z and $7 if n>0 and between —} 


m7 and 87 if n<0. 
Similarly for the argument of ce—im. When m > 0, n < 0, the result 
(16) means that either integral represents a function of ¢ analytic 


throughout the whole c-plane except possibly for c = 0. 
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In (16) put c = i/27(w.—w,), and replace x in the first and second 
integrals by 27a—22nw,, 27x—27ma, respectively. Hence 


is 4) 
e2mrixr dx 


(a—nw,)* (2 —Nevg) |’ 
. Wa— Wy 





nw, 


enix dx 


I(s) (x—meo,)!-* [“— 


Wa—W, 





Mr 


where, if n > 0, (2—nae) 
P— NW 
2) — Ia, 


3. v 
— jn <arg 


and if n< 0, , 
. U(x%—Nas) _ 5 
—37 < arg < 57, 


and similarly for arg(2—ma,). The paths of integration are lines 
from nw,, mw, to infinity parallel to the positive real axis of x and 


arg(x—nw,) = 0, arg(x—mw,) = 0 


on the respective paths. The only restriction on w,, w, in (17) is 
that w, ~ we, not I(w,) > I(we), since we can take c = 1/27(w,—ws). 
Hence (17), (17.1) hold if w, and w, are interchanged. 

I now use (17) to prove (6) wherein I(w,) > 0, I(w,) < 0, and 


Ax’+ Ban+Cn? _ (x—nw,)(x—nwg) 
vA y i(w.—w,) ; 


The path of integration in the left-hand side of (6) is deformed 
into the upper infinite semicircle, if m > 0, and into the lower, if 
m <0. The singularities of the integrand in the upper half-plane are 
given by «= nw, or nw, according as n > 0 or n < 0 and cuts are 
drawn from them parallel to the positive real axis of x. Hence four 
cases arise according to the signs of m, n. Denote the integrand of 
(6) by x(m, n, s), and take first the case m > 0,n > 0. Then 


io 2] 


f x(m, n, 8) da = (1—e?7*s) x(m, n, 8) dx, 


nw, 


since, when 2 describes a small circle around 2 = nw, in the clock- 
wise direction, arg(z—nw,) decreases by 27i, while the integral 
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around the small circle with centre at nw, tends to zero as R(s) < 
So x 0 
f x(n,m, 1—s) dx = (1—e-®7"s) f x(n, m, 1—s) dx, 
—2 Mwy 
since R(s) > 0. The arguments of the integrands at x =o are 
simply argi/(w,—wz,), i.e. they lie between +47, since arg(x—na)), 
arg(x—nw,) are both zero at x= 00. From (17) 

3 : wo 
J x(n, ”, 8) dat = e~FH88—-D (257) 28-1 Ta | x(n, m, 1—s) dx, 


ted | Mwy 


as in there , . 
= i(x—Nws) i(a—nNwa) 
arg 2’ | = —7-+arg| - = ae 
Way— Wy, W1—Ws 


for the left-hand side must lie between —3z and }z and, as we have 
remarked, —4n < argi/(w,—w,) < 47, while from a figure we see 


> ~ 


that arg(a—nw,) increases from 0 to arg(w,—w,) as x goes from oo 


to nw,. Hence 


oO 


yr “1—¢ ; 
J x(m, n, 8) dx = (27)*8-1 . ri | x(n,m, 1—s) da, (18) 
8 
and this is (6). 
Take next m < 0, n < 0; then 


| x(m, n, 8) dw = (l—e-?7"s) f x x(m, n, 8) dx, 


a NW 


ea) 4) 


x(n, m, 1—s) dx = (1—e?7"*) x(n, m, 1—s) dx. 


-—— MW, 


But (17) gives 
r(i—s) ¢ 


| x(m, n, 8) dx = er22-2)( 277) 28-1 re) | x(n, m, 1—s) da, 


7 
bet mw 


since there 
i(x—nw. “ i(a—nw. 
arg ( Ts mi+arg ( 2) , 
Wa— Wy, W1—We 


because the left-hand side must lie between —}z and 37. Hence 
(18) still holds. 
Take next m > 0, n <0, then 


| x(m, n, 8) dx = (1—e?7*s) f x(m, n, 8) dx 


n 


2 
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from the upper semicircle, while 


ao n 
f x(n, m, 1—s) da = (1—e?7"*) J x(n, m, 1—s) dx 
—2 Mg 
from the lower semicircle. For, in the case of the upper semicircle, 
arg(x—nw,) increases by 27 when 2 describes the small circle around 
nw», While in the other case arg(a—mw,) increases by 27. But from 
(17), on interchanging w,, wo, 
2) 
—8s 
@) ) | x(n, m, 1—s) dz, 
8 


MWe 


J xm, n, 8) dx = (27)*8—1 


NW 


rl 
r 


since we see from a figure that 


In < arg i(2—nw,) arg i(~a—mw),) <is 
W1— Wo, W1—Ws 
for the respective paths arg(~—nw,) = 0, arg(x—mw,) = 0, and these 
satisfy (17.1), if w,, w, are interchanged there. This gives (18) again. 
Similarly for m <0, n> 0. 
Now consider what must replace (6) when J(w,) >0, [(w ) > 0. 
We recall that VA denotes that value for which R(A/VA) < 0,and so 


Azx*+ Ban+-Cn? _ (x—nw,)(x—nwy) 
vA of 1(w,—Ws) 





> 





4) 
e2mmix dx 
I(m,n,s) = ah | x(m, n, 8) 
—@ 


say. Unless mn > 0, I(m,n,s)= 0 if we deform the path of integra- 
tion into the infinite upper semicircle when m > 0, and into the lower 
semicircle when m <0. It suffices to assume m>0, n>0. On 
drawing the cuts from nw, and nw, to 0 parallel to the positive 
real x-axis, and noticing that arg(z—nw,) decreases by 27 when x 
describes a small circle around nw, in the clockwise direction, we have 
Pa Pa 
I(m, n, 8) = (1—e?7"*) | x(m, n, 8) dx +-e?7's(1 —e?7"s) J x(m, n, 8) dx 


NW Ny 


= ¢,1,+¢,l5, say, (18) 
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where J,, J, denote the integrals. Hence 


x 
I(n,m, 1—s) = (1—e-*7") J x(n,m, 1—s) dx +- 

Ms 0 
+e-27!s(] —e-27!s) f x(n, m, 1—s) dx 

mw, 
= ¢)1,+c5]9, say. 
In both of these equations arg(v—nw,) = 0, arg(z—nw,) = 0 when 
r= ©. 
From (17), on interchanging w,, ws, 


(D0 9 (1—s) t 
— e—7i(2s—1)(97,)28—1 I 
I,=e (277) T(s) 1 


because in J,, 


xX—NwW, .(*x#—nw, 
arg] - = —7-+arg2({— ; 
~ Lt(@y—g) @1— We 


t— NW, 


and so as required in (17), arg ( 


lies between —3z and 47, 
W1—Ws a 
for, as 2 goes from 00 to nw,, the argument of the left-hand side 
varies from arg 1/i(w,—w,) to arg(—n/i) or at most from —z to 


DL 


Also from (17), by writing the denominator of (m,n, 8), as 


r—Nw, |* 
(w—nw,)*| = — |, 
1(w,;— 9) 


o,_1 1 (1—s 
we have I, (22)""* a : 8) I}, 
Ps) * 
' X—NwWy ' “ 
for as x goes from 00 to nw, arg . “_ goes from arg 1/i(w,—ws) 
. 1(wW1— Ws 
to argn/i and certainly lies between —37 and 37. Hence, if 


0 < Ris) <1, and I(w,) > I(w,), 
I(m,n, 8) —(27)8-1 


a — (e87is + 9 37is)(e7ls__e wie) T,. 


m 


9 vir 
e2mrix dx 
where I. => 


~ X—NWs, 8° 
(v—nw,)*| ; : 
1(w;— Ws) 
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This result, which is the analogue of (6), also holds for I(w,) = I(ws), 
if R(w,) < R(w,), as can be shown from (19) by considerations of 
continuity. 
Poisson’s summation-formula applied to f,(s) in (7) gives, as in 
R.S., when R(s) > 1, 


me _ = t e2mrix dx 
fils) = . 3 ze | (Aa?+ Ban+Cn?)*° 


Integration by parts shows that the series converges absolutely for 
R(s) > 0. The terms with mn < 0 vanish, and so 


since, in the series, n, m may be used instead of m, n as the variables 
of summation. Hence from (18), if 0< R(s) < 1, 
ae —8s) 


fils)—(277)*> y fill 


wn oa 


e2nriz dx 





wo 
- —(g37is_| g—87is)(e7is__ g—mis 2 4 
(c87i8 4 ¢-Sri8y( ) oc 
n=1 m=1 | (x—nw,)* | ————~ 
1(w—@g) 
MW, 


which, no longer assuming VA = 1, we can write as 


fy(s)—(2m VA) Tey fll—*) 


e2mnnir dx 
3zis Suis) wis___p—is 1-—2s 
— (e87'S-@ )(e eo") S S m | (de?+ Br+O)” 


n=1 m=1 


whee R(A/vA) <0 and I(w,) > I(w,), 0< R(s) <1. By expressing 
the integral in terms of the integral taken around the cut from w, to 
©, e.g. around the circle |z—w,|=e and cuts from w,+e to 0 
parallel to the real axis, the new form of the result is seen to be valid 
for bs values of s. 

. Another method of finding results such as (17) and (18) is 
sugge aan by using a differential equation. Thus 


e"' dt 
v= | iar (20) 
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with suitable limits for the integral, is a solution of the differential 


equation 
2(s—1)—= dy = 2 y+ dy ; 
dx dx? 


dy ite” dt L ms 1 
ae = ellxd . 
dz }(1—#)* 2(s—1) (1—#)41 


dy__—1 wel dt «x  f(1—é#)e™ dt 
(1—#)*-1-&(s—1)}  (1—#) 


if e(1—#?)!-* vanishes at the limits of integration, or returns to its 


dx  2(s—1) 


> 


initial value after describing the path of integration. And so, if the 
differentiation under the integral sign is legitimate, 


—1) 4 = (y+). (21) 
dx dx? 


But, if y=f(x,s) is a solution of this equation, another is 
y = x*8-1f(v, 1—s). For put y = 2-12, then 


lz om 
Me -28—1 ¢ +(%e y 228—2 
%e—1)| x 1 (2s—1)a | 


3 we + 2(28s—1)x*s- 1 + +(2s—1)(2s —2)2%-*2) ; 


da? 


x ( 1 


This is the same as (21) with s replaced by 1—s. Hence 
1-If(, 1—s) = Af (x, 8) + Bfy(a,8), 
where f(x,s), f,(x, 8) are two integrals of the type (20), and A, B are 
independent of x. The difficulty is now to find A, B, but there is no 
need to pursue the question. 
§ 5. The properties and functional equations of Epstein’s* zeta- 


* ‘Zur Theorie allgemeiner Zetafunktionen’, Math. Annalen, 56 (1903), 
615-44, and 63 (1907), 205-16. 
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functions can also be deduced by the methods previously explained. 
First, the summation-formula for n variables takes the form* 


® 
> F(x, Lg.+-%y) 


—-2 


Qmihrit... + 
e2riirit... itn) F(t, ty...t,) dt, ...dt,, 
_ oO 
where the summation refers to all integer values of x,,2%9,...,2,, 
provided (i) that 
Ortaat.-- +n W(x, Xe... 2, 
OxG).. . Ohm 

for a, @g,...,@, = 0, 1, 2 (excluding a,=a,=...==a,=2), are all con- 
tinuous functions of 2,, %,...,x, for all real values of 2, %,...,%,; 
and tend to zero when any variable tends to -+-00 regardless of the 
values of the other variables; 


(ii) that GP-S Fes 5 01 Bq) 
z a) D2 
Ox, Oz... Ox? 

Or F(x 


1 
2 2,2 
2 


Oxy Ox 


v — 
n) and similarly for 
n 


is an integral with respect to x, of “ 


the other variables; 


(iii) that 
a a 
|: : | > = aa a Ate Ly) dx 
2 OxF1 , . . 02% | nes 
—” x» Or =9,1 1 n 


converges. 
When n= 2, my proof shows that the right-hand side of (22), on 
integrating the general term by parts and writing 
20 A 
e2nrix e2nmix 
P(a 
2n7i)?’ 2nmi ’ 
= 
can be written as 


= 2 4) 


| | (Fe Per Gat PO ye tPOIPW 


ar 
_—.}dady. (22.2) 
”) 


CXL O° 


ey 


—-o —@ 
* *Poisson’s Summation-formula in several variables and some applications 
to the Theory of Numbers’, Proc. Cambridge Phil. Soc. 35 (1929), 412-20. 
695 H 
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Mr. Davenport points out to me that the identity of this with the 
left-hand side of (22) is seen more easily by using the formula 

Ns; | Ne Ne 

>’ F(n)= [ Fat — [ ROP dt, 

Ni N ie 
where NV, < N, and the dashes denote that, if NV, is an integer, the 
corresponding term is reckoned as }F(N) and similarly for N,. For 


Me Me 


Ms 
>” F(m,n) = | F(«,n) dx — | P, (a) 


m=M, 


OF (x, n) 


Ox 


dx. 
My WA 


Sum for n and note that the differentiation under the integral sign 


‘ 7 ‘ 7 a9 
A OF (x,n) . a : é oF 
arising from is legitimate, since ; are con- 
: Ox ' dx Oxdy 
tinuous. This gives 
Ng Ms 
" >” F(m,n) 


n=N, m=M, 


F 


Ma Ns 
-| | | Fe) - P,(y) : + P,(x)P,(y) be “ey | en 
ox 


¢ c 
OX oy oy 


My Ny 
Make M, > +, etc., and integrate by parts and we get (22.2). The 
proof is of course perfectly general. 

§ 6. Suppose now that 6=¢(x, x,...x,) is a homogeneous func- 
tion of dimensions A > 0 in 2, %,...,2%,, that 


(24, Le, « + +» Bq) = G(—2q, Ley. +3 —- Ve); 
that, for all real values of 21, 25, ...,%,, ¢ and all its derivatives exist 
and are continuous functions of x,, ...,2,, and that ¢ has a positive 
real part i.e. R(¢) >0. Finally, excluding 7,—2,=...=27,=0 and 
putting r2 — gt... +22, 
let |d|/r*, R(d)/r4 have positive lower bounds. 
Then a generalized zeta-function is defined by 


me (23) 


where the notation means that the summation is extended to all 
integer values of 1,, /,,...,1,, except simultaneous zeros. The series 
converges absolutely and uniformly for 


R(s) > n/A+6 > n/A, 
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since it may be compared with 
wo 
>’ G@+8...+12)-™. 
I=-—@ 


The summation-formula does not apply to the series (23), but it 
does to the series 


S [bly ly ...0,)He]l* = Fh ly... 1,) say, 


l=—@ 
where ¢ is a constant with R(c)>0, and R(s)>n/A. For we see 
that the conditions of (22) are satisfied by F(l, 1, ...1,) on noticing 


that, for Py = 0, 1, 2, 3..., o=1, 2,...,%, 


sp, | Mit FH a(ary Ly». + Ty) 
n 


pAt+pyt--- be 23.1 
dal... Gache ( 


is bounded for all real values of x,, 2, ... x, because of the con- 
tinuity and homogeneity of ¢ and its derivatives. 
Suppose finally that 


AMA [ ... f expl2ai(tyzy+...tt,t,)—mb (ty --- ty) dy... dt, 


= exp[—7@(x, 2... 2,)], (24) 


where ® is homogeneous of dimensions A with the same properties 
as $, Ay is independent of x, ..., #,, and A! “ denotes a specified root. 
This formula holds if ¢ and ® are interchanged* and Anat == i. 

I shall prove that ¢(s), defined by the series for R(s) > n/A, is an 
analytic function of s throughout the finite part of the s-plane except 
possibly for a simple pole at s = /A, and that it satisfies the func- 
tional equation 
n2s—nid Tg -+-n/d) 

AUA I(s) 


where (s), when R(s) > n/A, is defined by 


(s) = @(—s-+n/d) 


rs x, Xs 


* Cf. Muntz ‘Der Summensatz von Cauchy in beliebigen algebraischen 
Zahlkérpern und die Diskriminante derselben’. Math. Annalen, 90, 1923, 282. 


H2 








96 L. J. MORDELL 
is a quadratic form the real part of which is a positive definite form. 
Then A = 2 and 


2 


n 
Oe, ...%,) = > > 4&4 
1 


8? 


where A, A,, is the cofactor of a,, in the determinant 
Ag = |Q,4]; 
and Al? denotes the root with a positive real part. 
Let p, q be two positive constants and put 


2) 


¢,(s) = = [po(l, ...l,,)+qb(m, ...m,]-°, (25) 


lm=— 
where the argument of the expression in the brackets lies between 
+147, and the summation extends to all integer values of /,, ..., l,, 
M,,...,m, except 1,=...=l,—m,=...m,=0. It will be seen that 
the series converges absolutely if R(s)>2n/A. The terms with 
m,=...=m,=0 give p-*d(s). For any other given set m,, ..., m,, 
apply the formula (22) to the summation in /. Then 


2) 


ah 0 © { e2miitit...+ln) dt, ... dt, 
$1(8) = p*9(8) + P3 > | an | sen ...t,)+¢qb(m,...m,)}* 


Isolate the term with /,=/,...—1,—0 and write in it 
[q¢(m, ... m,)/p]}*t, 
for t,, ete. Since R(s) > 2n/A, 
$1(8) = p-*h(s) pgs" b(s—n/A)K ($, 8) + 
e Ss" = 2, | e2mi(hti+... +Hntn) dt, a dt, , (26) 
i. afta t. [pblt ...t,)+qp(m, ...m,}° 


wherein the combinations /,=/,=...—1,=0 and m,=...=m,=0 
are both omitted from the summation. Also 


Dn 


lt, 


: 7 (* dt, ..s4 
Kg) =[...| oy new t,)+1]}*’ 


the path of integration for each variable being the line joining 


+d(m,...m,,)"A 00. 
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Denote the last integral in (26) by I(l,m,s,¢). Clearly 


io 3) 4) 
a) 


| expl2nitits+. .-+1,t,)— 


0 


(s) 
—ntpd(t, ...t,,)—mtgd(m, ... m,,)]t®1 dt, ... dt, dt. (27) 


Replace t, by t,(tp)-, &c., where the roots are positive numbers 


(I, m,8,$) = [--. 


—@ 


and apply (24), whence 


2 


8 m—nlA 
I(l,m,n, 6) = re) ‘AYA | expl—matg(on ...M,)— 
0 
—7®(l, ...1,,)/tpys-1-" dt, (28) 
the justification of the interchanges in the order of integration being 
quite easy. 

From (27), on taking /,=...=1l,=0, p=q=¢(m, ... m,)=1, and 
provided that the path of integration fort,,¢, ... can be deformed* into 
that in K(d, s), we have 


fo) 


K( ,8) == ar eee A dt 
¢ AYsT(s) 


0 
—_ wMAT(s—n/A) 
* AY I(s) 
Next, (28) shows that 
I(l,m,n,¢) = Ofexp[—M.RO(I, ...1,,).Rb(m, ...m,)]} 
where M > 0 and the constants in O are independent of 1,, ..., l,, 
M,,..+, m,, the result holding uniformly in s for bounded s. 
For all t, js—1-n —_ Ofexp(5t+8/t)], 
where § > 0 is any fixed arbitrary constant smaller, say, than both 
aq R[d(m, ...m,,)], aR[O(l, ...1,)]/p. 


This is permissible since R[¢(m, ... m,,)] has a positive lower bound 
for integer values of m,, ..., m, and similarly for O(/, ... l,). The 
inequality follows, since 

ft+g/t > 2v(f9) 


for f>0,g9>0,t>0, 
and | e-vt+10 dt converges, if y > 0. 
0 


* As is the case when ¢ is a quadratic form. 
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This shows that the series in (26), whose convergence is as yet 
undemonstrated, converges absolutely for all values of s, and defines 
an integral function of s. This result is also easily seen by writing 
2ril, e? 74h — d(e?7hh), ete., when /, 40, and repeatedly integrating 
by parts with respect to the variables. 

The left-hand side of (26) is a symmetrical function of p and q. 
Interchange p, g and subtract, and we have 


(p-*—q-*)$(8) + 
+ (p-ng-*+"A—q-" psn) $(s—n/A) K(f, 8) = Gy(s), (30) 


where G(s) is an integral function of s. This is really a recurrence- 
formula for ¢(s) from which we deduce that 


(p- n Ag-* +n A_q-" \p-s +n \\b(s—n, A)K(¢, 8) 
is an analytic function of s for R(s) > n/A. 


Change s into s+-n/A and note (29); then 


(p-"9q-*—q-™\p-*) (5) F(s)/P(s-+n/) 
is regular for R(s) > 0. Since the first factor has a simple zero when 
s=n/X, we see that ¢(s) is regular for R(s) > 0 except possibly for 
a simple pole at s=n/A. The functional equation (31) then shows 
that ¢(s) is regular for R(s) < 0, no difficulty arising at s = 0 since 
®(—s-+n/d)/T'(s) is reguiar at s = 0. 
The question whether s=7/A is a pole can be often answered 


by Landau’s theorem on the singularities of Dirichlet’s series with 


positive terms. 
Now write 1/tpq for ¢ in (28) which becomes 


ins) 


n\ 
. | expl—nat\ .-b,)— 


—nd(m, ... m,,)/tp}t-s-1*" dt. 


an 
I (l,m, 8, $) = P(s)A's 


Substitute —s-+2n/A for s and ® for ¢, and we have 


AMAT (s) AWA T(—s-+ 2n/d) 


A nid ps 2njA ge—nia 


I(l,m,8,¢) = I(m,l, —s+2n/A,®) 


7 p-" m8 +2 


or 

I(l,m, 8, ¢) 

, ALA T'(—s-+ 2n/d) 
Als I(s) 


28—2n (pq)-8*"AI(m, 1, —s+2n/d,®). 


= 7 
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Then (26) gives at once 
$1(8)—p~*h(s)—p-"q-##" K(f, 8) b(s—n/A) 
ALA T(—s+2n/aA 
eT a [Oa —04 2a) 
— p*-2"AD(— s+ 2n/A)— pgs "AK (®, —s+ 2n/A)O(—s+-n/d)] 
or say ¥,~¥ 0 pat V.—V,—V,. 
Note that V,—V, is a symmetric function of p and gq, while 
p'(Vo—Ve), pgs "*(V,—V;) are each the sum of a term independent 
of p and q together with a power of pg. By interchanging p, g, we 
see that each of these three terms V,—V,, &c., is zero, since 


(p—q*)/(prrgs-°—qr ps") 

is not a function of pg for general s. Hence 

», , AYA I'(—s+2n/a ae ; 
ra Fay eatery 8-420), 
2t—2nid AVA T'(—s8+2n/d) 

A's I'(s) 
From (30), this can be written as 

a 

¢ (s) 

For the quadratic form, this is Epstein’s* well-known result. The 
same method gives his corresponding formulae when the general term 
of the zeta-function has an additional exponential factor of the usual 
type. 

§ 7. There still remains the study of the zeta-function arising from 
quadratic forms of at least three variables, which, when multiplied 
by any constant factor, have indefinite forms for their real parts. 
This depends upon considerations of an entirely different character 
from those which have preceded. For now the zeta-functions either 
do not exist or the question of their existence is bound up with the 
arithmetical nature of the coefficients of the form. But this involves 
a field wherein not much is known. Take, for example, the case of 
four variables in (20) when 

f = (ly ly ls lg) = py (ly Ly lg lg) + tyho(l, Lys 1), 
—— hy = Al} +agl3+asl3+-a,lj, 
the = bl} +byl3 +6513 +b,i, 
* loc. cit., Math. Annalen, 63 (1907), 207. 


2s— 
— 75 2n/ 


(8) 


d(s) = K(®, —s+-2n/A)®(—s-+n/d). 


@(—s-+n/d). (31) 


57°805A 
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and a, @, ..., 5,, bs, ... are rational integers. It seems that they 
can be chosen so that ~, = %, = 0 for an infinity of integer values of 
l,,...,l,. For eliminating 1,, 
A, l?+-A,13+-A, 13 = 0, say. 
The integer solutions of this equation are given by /,=1,=/,=0, 
and also, when certain conditions are satisfied, by 


1 = xi(P 9) lp = x2(P, 4); ls = x3(P. 9); 


where x;, x2, X3 are homogeneous binary quadratic forms finite in 
number, and p, g are any integers. Then 

9 

i= x(P-9), 
where x is one of a finite number of homogeneous binary quartics. 
But for many such quartics, there are an infinity of solutions in 
integers /,, p, g, and in this case the zeta-function cannot exist. More 
generally it may be possible that ¢ can take the same value for an 
infinity of integer values of 1,, /,, l,, 1,. Similarly, for more than four 


=...=0. 


> 


variables, since we need only take /; 

For the case of three variables, take the special case of (20) when 

d = d(l,m, n) = m?—2l?+-1(n?—3l?). 

First of all the application of Poisson’s formula is not valid, since 
d(x, y,z) becomes zero when y = xv2, z= 2v3. But apart from this, 
it is not easy to ascertain whether the series defining the zeta-function 
converges or not. 

Thus, for given integers M, N, there are at most only a finite 

} > . 
number, say 7'(M, N), of integer solutions of 
m*?— 2]? = M, n?—3l? = N. 
For, if p= mn, : , ee 
I p? = (2/?-+ M)(3l?+-N), 

and such equations have only a finite number of integer solutions 
in p, 1. But as very little is known about 7(M,N), it is not very 
useful to compare (20) with 


> >’ TUM, NP +N) Ane, 


— wr 


Take, for another case, 


dh = (aym+b,l)(agm+byl)+i(e,n+d,1)(con-+d,l), 
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where the constants a,, 6,, ... are irrational numbers. Apparently 


they can be selected so that 


lagm+b,l| < |l|-}, |con+d,l| < |l|-* 


for an infinity of integer values of 1, m,n. For these values of 1, m, n, 
¢ is clearly bounded, and so > ¢-* diverges and the zeta-function 


does not exist. 
These illustrations are sufficient to show that a new type of 
question has arisen. 





A NOTE ON SERIES OF SINES 
By A. ZYGMUND 
(Received 27 March 1930) 
§1 


1. In this note we consider some properties of the trigonometric 
sine-series os 
> 6, sin nx (1) 
n=1 
in the interval (0, z). It is well known that the functions sinz, 
sin 2x,..., form a complete set of orthogonal functions in this interval ; 
and therefore every function f(7) integrable in the interval (in the 
sense of Lebesgue) can be developed in a series of the form (1), with 
o% 
= | f(t) sin nt dt. (2) 
0 
These formulae may, however, preserve their sense even when /(#) is 
not integrable. It is sufficient to suppose that f(¢)sin¢ has this pro- 
perty or, in other words, that ¢f(t) and (7—1t)f(t) are integrable over 
the intervals (0, $7), (47, 7) respectively. ‘In this case the coeffi- 
cients b,, do not necessarily tend to 0, as is shown by the example of 
the function f(7) = } cot $x, for which b, =b,=...=1. It is, how- 
ever, easily proved by partial integration that in the general case 
b, =0(n). 
When /(/) possesses the properties just stated and b,, is defined by 
(2), we shall say that the series (1) is a Fourier sine-series. 
2. Series of the form (1) do not possess some properties of general 
trigonometrical series 


we 
> (a, cos nx + b, sin nx) (3) 
v=1 


39+ 
7 


considered in the interval (—7z, 7). In regard to the series (3) we 
have, for example, the following theorem: + 
(SB.) If a series (3) converges in every point of the interval (—z, 7) 
to a (finite) sum f(x) satisfying an inequality f(x) > (x), where d(x) is 
integrable, then f(x) is itself integrable and (3) is the Fourier series of 
t+ Due to H. Steinhaus and 8S. Banach. Cf. H. Steinhaus, ‘Some properties 


of trigonometric series (Polish)’, Cracow Rozprawy, 56 (1916); A. Zygmund, 
Math. Zeitschr. 25 (1926), p. 278. 
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/(x). The proposition is still true when the series (3) fails to converge 
in an (at most) enumerable set of points E.+ 
An analogous theorem holds for the cosine-series 


ao 
lag + > a, cosnx 
=1 


n 

considered in the interval (0, z). It is a trivial deduction from (SB). 
But the corresponding theorem for sine-series is false. In fact, the 
=. sin nx 
z log n 
n=2 
converges everywhere and its sum is bounded below without being 
integrable. We can, however, prove the following theorem. 

(A.) If the series (1) converges throughout the interval (0, 7), except 
possibly in a set E, to a finite sum f(x), and if, except in E, 

fx) > $(a)/sin a, 

where ¢ is integrable over (0, 7), then f(x) sinx is integrable and (1) is 
the Fourier sine-series of f(x). 

3. We prove our theorem by reducing the sine-series to a cosine- 
series. Multiplying the N-th partial sum of the series (1) by sinz, 


series 


we get 
N N 
sin x > b,, sin na = $¢9+ 2m cos nx — by cos(N+1)a—by,, cos Na, (4) 
1 


n n= 


where Co = 6,, Cn = $(On 44 —Oy-4); (n= 1,3,...0). . @ 
Now, since b,, tends to zero with 1/n, the series 
d¢o+ 2, on COS NX (6) 
behaves in the same manner as the series (1) multiplied by sinz. 
More exactly, the difference of these two series converges uniformly 
to 0 throughout the interval (0, 7). Therefore the series (6) converges 
to f(x) sin x wherever (1) converges. It follows from (SB) that f(x) sin x 
is integrable and that 
2 7 2 7 
Se = | f()sint dt, C, == | f(d)sintcosnt a (n=1,3,...) 
~ ) 
4 


Putting, then, bF = 


[fesin nt dt, 
0 


7 


+ We shall throughout the paper denote by £ a finite or enumerable set 
of points. 
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we obtain from (7) and (8) 
Co = OF, Cy = (Oh 44-574) 

which, with (5), gives us =o 

4. The conclusion of the theorem (SB) is true under more general 
conditions. It has been proved,} for example, that when the series 
(3) is perhaps not convergent, but its upper and lower sums f and f 
are finite (except, possibly, in a set F), and f exceeds an integrable 
function ¢(x), then (3) is the Fourier series of f. It has been shown 
also that this generalization remains true when we define fand fas 
Poisson’s upper and lower sums, i.e. by ad 


f(x) = lim [aot 


00 
~ 
" r—>l1 - 


(a, cos nx + 6, sin nar : 


1 


n” 

a ee D 

f(x) = lim [$a0++ > (a, cos nx + 6, sin nx)r"|, 
r—1 n=1 


provided that a,, 6, tend to zero.t Using the same method of proof 


nm 


as above, we obtain the theorem 
(A,.) If b,, + 0, if (except, perhaps, in E,) the functions 


a) mi 0 
f=lim > b,sinnzr", f=lim > b,sinnar” (9) 
= ro>ji"=i1 r—>1n 1 
are finite, and f(x) > 4(«)/sinx, where ¢ is integrable, then f(x)sinx is 
integrable and (1) is the Fourier sine-series of f. 
5. (B.) If b, + 0 and 
Ab, — b,—b,, +1 Sm 0 ’ (10) 
then the series (1), which converges everywhere, is the Fourier sine-series 
of its sum. 
In fact, writing 
nt 
Q ») . 13 . 
8,(x)= > b,sinra, 
n=1 , , 
: F ‘ cos 5% — cos(n-+-3)a 
t, (x) = sina + sin 2a+...+ sinna = : : (n-+3)2 , 
2 sin $a 
we have, for 0<2< 7, 
N+k Nik 
b,, sin Nz = 3 t,,(a)Ab,, —ty by 1 wttner bys k+1> 
n=N+1 n=N+1 


1 aah 
N+ k—58N S p 2 1 Ab,, a i by utby, k+1} = 


~ sin $x \n=N3 


Sy tk Sn = 
9 

2b 44 4bn 44 
sinjx ~ sinx 


+ By Steinhaus and Banach (loc. cit.). 
t A. Rajchman, Prace Mat. Fiz. 1919; cf. also A. Zygmund, loc. cit. 
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From this last inequality it follows that (1) is uniformly convergent 
in the interval (e, 7), for every positive number e, while 


~ 
> 6,, sin nx sin x 
n=1 
is uniformly convergent in the whole interval (0, z). Accordingly, 
denoting by f(x) the sum of the series (1), we see that the series (6), 
whose coefficients are given by (5), is the Fourier series of f(a) sin 2. 
Thus, reasoning as in the case (A), we obtain the desired result. It 
is evident that in this theorem the condition (10) may be replaced 
by > |Ab,, | < @, 


§ 2 
6. The conditions imposed on b, in the theorem (A,) may be 
weakened, provided that certain restrictions are made as to the 
behaviour of (1) in the exceptional set Z. 
(A,.) Suppose that the series (1) satisfies the conditions of (A,), except 
that b, = O(n*), where k <1. Then if 


lim (1—r) S b,, sin nav, r” = 0 (11) 
r—>1l1 n=1 
in every point x» of E, the conclusion of (A,) remains valid. 

7. The proof follows the same lines as that of (A,). Let us suppose 
first that the remainders by, cos Nx, by cos(N-+-1)x in the formula 
(4) tend to zero by Poisson’s method in the open interval 0 < 2 < z, 
i.e. that 


lim (1—r) 5 £,,(x)r" = 0 (0<2 (12) 
n=1 


r—>1 


where &,(x) denotes either 5, cos(n+-1l)a or b,,,,cosna. Then the 
Poisson’s sums of the series (6) are fsina and fsinz, and are finite 
for the points belonging to the set Z. At any point x, belonging to 
E we have in any case 


2) 
lim (l—r)(4¢9+ > ¢,, cos nay r") = 0. (13) 
r—>1 n=1 


Taking into account the form of the coefficients c, in (5), and the 
fact that b, = O(n*)=o(n), we find that the relation (13) holds 
for %=0 and a,=7z. It follows that (6) is the Fourier series 


+ A. Zygmund, loc. cit. 
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of f sin x. Accordingly, as in (A,), the series (1) is the Fourier sine- 
series of f. 

8. It remains now to prove that the relation (12) results from our 
hypotheses. To this end we observe that (12) holds, if 


lim (1—r) > 6, sinna r” = 0, (13a) 


7 1 n=1 


lim (1—r) > b,, cos nx r" = 0, (0<a<7) (13d) 
r—1 n=1 


hold. (13a) of course follows from (11) and the finiteness of f and f 


outside H. To prove (13b) we may reason as follows. From the 
hypotheses of the theorem (A,) it results? that the series (1) may be 
integrated term by term in any interval (e, 7—e«), and the sum so 
obtained converges to the integral of f(x) [or f(x)]. From the con- 
vergence of the series > (b,,/n)cosnx it results that the sequence 
b, cosnx (n= 1, % ) is summable (C, 1), and also by Poisson’s 
method, to zero, and the relation (130) is proved. 


9. (B,.) If 6, =o(n), > n|A%b,,| < 0, (14) 


then the series (1) is uniformly summable (C, 1) in every interval (e, 7) 
to a function f(x) such that f(~) sin x ts integrable. Moreover, the series 
(1) is the Fourier sine-series of f(x). 

We observe first that, because of our hypotheses, qg,, = Ab, = 0(1). 
We consider as above the difference of the series (1), multiplied by 
sin x, and (6). This difference is summable (C, 1) to zero if b,, ,, cos nx 
and b,, cos(n+1)a tend (C, 1) to zero. To show that this is true, it is 
sufficient to prove that b, cosnxz and 6, sinnx tend (C, 1) to zero, 
that is, that w 

>= b, e"* = o(N) (O<a< 7). (15) 

n=1 
But this is an immediate consequence of Abel’s transformation and 
the fact that b,, = o(n) and Ab, > 0. Moreover, it is evident that (15) 
holds uniformly in any interval (c, 7). Taking into account now the 
fact that 

Q= 3 (On 41—5,,-1) emg. 3(Qn-1+9n)s 
we see that the coefficients of the series (6) satisfy the conditions 

¢, = 0(1), > n| A%c,,| < 00. 


+ A. Zygmund, Math. Zeitschrift, 24 (1925), p. 93. 
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Hence, by a well-known theorem of W. H. Young,} (6) converges 
uniformly in any interval (e, 7) to a function F(x) and is the Fourier 
series of F. The series (1) is then uniformly summable (C, 1) in the 
interval (e, 7) to a function f(x) and evidently F(x) =f(x)sinz. We 
then have (7), and using (8), we obtain by the same argument as in 
(B) that b, = b* 


n° 


t W. H. Young, Proc. London Math. Soc. (2) 12 (1913), pp. 41-70. 











A PROBLEM CONCERNING THE HYPER- 
GEOMETRIC EQUATION 
By E. G. C. POOLE 
[Received 25th April 1930} 

1.1. Certain well-known homogeneous linear differential equations 
are most easily solved by utilizing some simple algebraic relation 
between the solutions, known in advance. Thus Lamé’s equation of 
integral order has in general no periodic solution; but it has two 


quasi-periodic solutions whose product is an elliptic function. The 


latter is not difficult to calculate in any given case and the equation 
can be solved by a quadrature. 

1.2. Schwarz’s investigationt of the cases of complete algebraic 
solubility of the hypergeometric equation is too celebrated to require 
detailed comment. There we have a group of solutions, every sym- 
metric function of which is a rational function of the independent 
variable ; we are tempted to ask what would happen if only one such 
symmetric function need be rational. Several investigations on these 
lines have been made. Markoff{ solved various problems, in which 
he assumed the validity of some simple linear or quadratic relation 
between two solutions of the hypergeometric equation, or between 
a solution and its derivative. Among these was the problem of deter- 
mining the conditions that the product of two solutions should be 
rational. His method was algebraic, and consisted in forming an 
equation of the third order satisfied by the product of two solutions 
of the original equation, and in examining when this equation has 
a rational solution. 

1.3. A more direct method was devised by van Vleck§, who 
extended it to the general Fuchsian equation of the second order, 
with any number of singularities. In the present paper I shall 
restrict myself to the equation with three singularities (Riemann’s 
hypergeometric equation), and I shall begin to examine cases where 
the product of three or more solutions is rational. If n is the number 
of solutions, the case n = 3 is easily treated on the lines of the known 

+ Crelle, 75 (1873), 292-335; reprinted Ges. Werke, ii. 211—59. 

t Math. Annalen, 28 (1887), 586-93; 29 (1887), 247-58; 40 (1892), 313-16. 
Comptes Rendus, 114 (1892), 54-5. 

§ Am. Journal of Math. 21 (1899), 126-67. 
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cases m=1 and n=2. There are a number of easy lemmas 
applicable to all values of n, but the detailed results for n = 4 lead 
to rather heavy algebra and are still incomplete, so that they will 
not be included in the present paper. 

2.1. We begin by a trivial simplification. In the equation of the 
second order, the most general solution is a linear combination of 
two independent solutions. If we consider a set of m solutions and 
a homogeneous polynomial of degree n in those solutions whose 
coefficients are constants, the polynomial can be replaced by a 
quantic of order n in the two fundamental solutions. On factorizing 
the quantic, the original function of m solutions becomes a product of 
n solutions. Results obtained for such a product admit of generaliza- 
tion for a polynomial. 

Let Y,, Y., ..., Y, be solutions of Riemann’s hypergeometric 
equation, represented by the symbol 


A BC 
Ply PB, wy XY], (2.11) 
“2 Bz Ye 
where the six exponents must satisfy Fuchs’s relation 
Oy +a+f8i+Potnut7e= 1. (2.12) 


ProsBLemM. Under what conditions is the product Y, Y.,...Yy 
a uniform function of X ? 
We make the classical change of variables 


___ (X—A)(B—C) _ (X—A\-*(X—C\ ny 213 
* = (X—ByA—C)’ oo xX—B : 


and obtain the symbol 


0 00 1 
r( 0 (By+4+71) 0 vw], (2.14) 
(%,—o4) (Baty +71) (Y2o—11) 
which is that of an ordinary hypergeometric equation, satisfied by 
F(a, b;c; x), where 
a=a,+8,+7, b=a,+h.+y1, c=1+a,—ag. (2.15) 
We may therefore work with the usual Gaussian form of the 
equation 
2(1—2) 44 {e—(a-+b+ 2 aby = 0, (2.16) 


and the question is modified as follows: 
3695 I 
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ProsieM. If ¥;, Yo, ---,Y, are solutions of (2.6), under what con- 
ditions is the product y, Y2...Y, Of the form 

xP(1—a)* x (a uniform function of x) ? 

2.2. It is convenient to have a short name for an analytic function 
whose logarithmic derivative is uniform, so that two branches ob- 
tained by analytic continuation from one another only differ by a 
constant multiplier. I shall venture to call such a function mono- 
dromic, a term now rarely employed, but formerly meaning a function 
which is itself uniform. We therefore ask the question, when has 
the hypergeometric equation a set of m solutions, whose product is 
monodromic ? 

2.3. Lemma. Let (y;), (y%) be two sets of n solutions of a homo- 
geneous linear differential equation of the second order. If we have 
Y1 Yo +++ Yn = Yi ys ... YX, then each set is a rearrangement of mul- 
tiples of the other. 

To prove this, we express each side as a binary quantic, in two 
independent solutions u, v. If the two quantics have the same coeffi- 
cients, they admit the same linear factors, and the theorem is proved. 
But if the coefficients are different, we can form a quantic in u, v, 
which vanishes identically, although its coefficients are not all zero. 
This cannot be true, if w, v are linearly independent. 

2.4. The foregoing lemma enables us to reduce the problem to the 
case, where any pair of the solutions is linearly independent, and 
each solution occurs to the first power only. For, if the conditions 
contemplated in the problem hold, we may group together the »,. 
solutions occurring to the power k. On examining the effect of 
analytic continuation about a circuit, each solution must become one 
of the set occurring to the same power as itself. Thus each separate 
set is separately monodromic. Taking the k-th root, and observing 
that the product is free from singularities except at x = 0, x = 1, and 
x = 00, we see that the product of the solutions of the set taken to 
the first power is also necessarily monodromic. 

3.1. The solution for n = 1 is given by both Schwarz and Markoff. 
The monodromic solution is one of the fundamental solutions at each 


singularity. By a trivial transformation, its exponent at x= 0 and 
at| «= 1 may be taken as zero. The solution is then regular and 
uniform throughout the finite part of the plane; and since 2 = © is 
a regular Fuchsian singularity, it can only be a pole. Thus the typical 
case, to which all others can be reduced, is that of the hypergeometric 
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polynomial F(a, b; c; x), where either a or b is zero or a negative 
integer. 

3.2. Can there be a second monodromic solution? In order that 
the second solution of (2.6) may be free from logarithms, it is neces- 
sary that the exponent differences 1—c and c—a—b should not be 
zero; special subsidiary conditions are required if either of these is 
an integer. The second monodromic solution, if it exists, must be 

Yo = x-“(1—x)°-4*  F(1—a, 1—b ; 2—c;; 2). (3.21) 
3y 3.1, the condition for the first solution may be taken as a equal 
to zero or a negative integer. The hypergeometric function in (3.21) 
is to be regular in the finite part of the plane and therefore a poly- 
nomial, which requires that b shall be a positive integer. The con- 
ditions are therefore that 
ax<0<1<8), [a, b integers] (3.22) 
and this makes x = 00 an apparent singularity. 

3.3. This solution fails if ¢ is an integer and if 2<c <b. Other 
integral values of c, however, yield a uniform solution (3.21) and the 
most general solution is then uniform. As the linear transformations 
of the hypergeometric equation are well known, we shall not attempt 
to enumerate all possible sets of conditions, as they can be reduced 
to the present case. A typical example is 


b<c<a<0, [a, b, c, integers] (3.31) 
where the solutions involve two polynomials 
y, = F(a, b, c, x), Y= v1“ F(a—c+1, b—c+1; 2—c; 2). (3.32) 


Solutions of this kind yield an infinity of trivial solutions of the 
general problem, since any polynomial in y, and y, is a uniform 


function. 
4.1. In sketching the solution for n = 2, we follow the methods of 
van Vleck rather than Markoff, in view of future developments. 
Suppose, then, that y, y. is monodromic, and consider the effect of 
analytic continuation about any closed circuit. Neglecting constant 
multipliers, y, and y, must either be interchanged or reproduced. If 
circuits about two singularities interchange y, and y,, a combined 
circuit about both of them (which is deformable into a circuit about 
the third singularity) will reproduce them. Thus one point always 
reproduces y, and y,, and, by a linear transformation, this may be 
taken as x = 00. It is obvious that if two points reproduce y, and Ye, 
12 
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then the third also reproduces them. Thus the singularities at x = 0 
and x= 1 will either both interchange, or both reproduce y, and yp. 
If all three singularities reproduce y, and y,, then y, and y, are 
separately monodromic, and we have the case discussed in 3.2. 
4.2. Suppose, therefore, that simple circuits about x = 0 or about 
x = 1 interchange y, and y,; then a double circuit about either point 
reproduces them. The points are fundamentally quadratic irra- 
tionalities, where the exponent difference is half an odd integer. We 
can adjust the exponents so that, instead of (2.16), we have the 
Riemann equation corresponding to 


( 0 00 ] 
Ptk(m+1) uk §(n+1) 2}, [m, n integers] (4.21) 
—km —hik —43n 

and on putting x = sin*0, we have the equation 

d*y 

doz 
an equation solved by Elliott and by Darboux.t We know two 
solutions y, = e*R(cot@) and .y,—e-**R(—cot@), where R is a 
rational function. These are distinct, unless k = 0, and 

Y1 Ye = R(cot 6) R(—cot @) 


is a rational function of 2 = sin*6. Some comments on this equation 


[m(m-+1) cosec?@+-n(n-+ 1) sec?6-+-k?]y, (4.22) 


by the present writer appear elsewhere.§ 

5.1. The characteristic property of a regular singularity is the ex- 
ponent difference; if this is a rational fraction, whose denominator 
(when expressed in the lowest terms) is m, then we may say that the 
singularity is of order m; any solution whatever of the equation 
becomes a multiple of itself after describing m circuits about the 
singularity, provided that m is greater than unity. If m is unity, 
the point is usually a logarithmic singularity, but may be an apparent 
singularity, in which case we shall say the order is 1. 

5.2. Lemma. Under the conditions of the problem, if n is greater 
than 2, each exponent difference must be real and rational. 

Using Riemann’s form (2.11), let P,,, P,, denote the two funda- 
mental solutions at a singularity. Since n > 2, and since each solution 
occurs only to the first order, the product y, y, ... y,, must contain 
at least one factor y,=AP,+BP,,, where A#0, BAO. After 

+ Elliot, Acta Mathematica, 2 (1883), 233-60. 
t Darboux, Théorie Générale des Surfaces, ii. 207-13. 
§ Journal London Math. Soc. (shortly). 
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s circuits, this becomes yf = Ae?" P, + Be? P, , and, for all 
integral values of s, this must be a multiple of one of the factors 
Yis Yor «++ Yn. Thus e?7(e2-™)s can take only a finite number of values, 
and accordingly (a,—«,) is rational. 

5.3. Corottary. The order of a singularity is not greater than n. 
For if (a,—a,) is a fraction in its lowest terms, with denominator m, 
the general solution has m distinct branches, which interchange 
around the singularity, each of which must occur among the n factors 
of the monodromic product. 

5.4. Lemma. Under the conditions of the problem, if one singularity 
is apparent (i.e. of order unity and free from logarithms), the other two 
must be of the same order m. 

Every solution is monodromic, with respect to a circuit about the 
apparent singularity, which is deformable into a sequence of two 
circuits about the other singularities. If one of these is of order m, 
consider a group of m solutions which are cyclically permuted by a 
circuit about this singularity; since the combined circuit leaves each 
solution unchanged (neglecting a constant multiplier), the circuit 
about the second singularity must reverse the permutation effected 
by a circuit about the first, and this requires that the cycles should 
be of the same order m. 

6.1. The case n = 3 can now be rapidly disposed of, by enumerating 
possible combinations of singularities, none of which may be of order 
greater than 3. 

6.2. THEoREM. There can be two or more singularities of order 3, 
only when the conditions of (3.2) are satisfied, so that two solutions are 
monodromic. 

Without loss of generality, let « = 0 and x = 1 be points of order 3. 
The monodromic product is evidently expressible in two forms, which 
become (on absorbing trivial coefficients) 

P+ ert 


a 


where P, = pP+4Fh,.. P _= rP +sP, ps—qr £0. 


a, a? 
Hence ptr = 1, g@?+s3 = 1, 
pqtrs=0, pq?+rs?=0. 
We have, since ps—gqr £0, 
pq = 9, rs = 0, 


and so p=s=0, or q=r=2, 
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which only differ by a change of notation. Taking the second form, 


we are left with p=1, g=1, 


and by a proper choice of notation, we may take 


p = iF q = & 
Thus P =P, and P= * 


Oy vv Og 
and we have two solutions monodromic at each of the two points 
x = 0, x = 1, and therefore also at x = 0. 
6.3. The above theorem and the lemma of 5.4 reduce the cases 
left to be examined to four, the orders of whose singularities are 


respectively (1,1, 1), (1, 2, 2), (2, 2, 2), (2, 2, 3), 


which we shall examine in turn. 

6.4. The combination (1, 1, 1) is that of the trivial case disposed 
of in 3.3, where every solution is uniform. 

6.5. Each of the three remaining cases has two singularities of 
order 2, and therefore admits of the reduction of 4.2, since the two 
points in question may be taken as x = 0 and x= 1. The remaining 
singularity may be of order 1, 2, or 3. The fundamental solutions 
at x = 00, after the transformation x = sin?@, become 

P;, = e R(cot 8), P;, = e~*® R(—cot 8), (6.51) 
where the original exponent difference tk is rational, with denomina- 
tor 1, 2, or 3. We have B,—f, = +1k. 

6.6. If ik is an integer, we have the case (1, 2, 2). The product 
Y1 Y2 Yz is monodromic at x= oo, if it is monodromic when @ is 
increased by 7, which is true for any homogeneous quantic in F,. 
and F;. The product is monodromic at «= 0 if it corresponds to 
either an even or an odd function of @. The conditions at x = 1 will 
be satisfied automatically. We can construct the following com- 
binations 

(P3+P3), Ps Pa,(Ps, Ps), (6.61) 
which show that no further conditions are necessary. 

6.7. If ik is half an odd integer, we have the case (2, 2, 2). If the 
product exists, it can be written in the form 

V1 Y2Y3= AP3.+ BP} P,.+ CP,,P3,+ DP%,,. (6.71) 

On increasing @ by z, we find that this must be a multiple of 

Liyt y3 y} = AP3 —BP3 P3,t+ CPs P3,—DP3,, (6.72) 
hence either A = C= 0, or B= D= 0. 
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As these only differ in notation, let us take B= D=0. We are 
left with the requirement that AP +CP;,P}, should be either an 
even or an odd function of 0, which is impossible. Thus the case 
(2, 2, 2) is excluded. 
6.8. If the denominator of ik is 3, we have the case (2, 2, 3). The 
product must be of the form 


"1 ¥2Y¥3 = AP} + BP3,, (6.81) 
in order to be monodromic at «=0oo. To make this an even or an 
odd function of 0, we require A = + B, so that the following com- 
binations satisfy all the conditions 


Yi Y243 = Pit P%3.. (6.82) 
The enumeration of the possibilities for n = 3 has now been com- 
pleted. 








ON THE EXPANSION OF THE EARTH’S GRAVITY 
IN POWERS OF THE SQUARE OF THE SINE OF 
THE LATITUDE 
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In the case in which the Earth’s surface is a surface of equilibrium 


as well as an oblate ellipsoid, Stokes has given an expression for the 
external potential correct to the first order of approximation.* But 
the process may be continued, as I have shown in two recent papers, 
where I have established the rigorous form of the expansion of the 
potential above mentioned and consequently the expansion of the 
Earth’s gravity in a series of Legendre’s functions. 

My object in the present communication is to give the expansion 
of the gravity in powers of the square of the sine of the latitude, 
such an expansion being the most interesting for practical applica- 
tions. 

Moreover, it is to be observed that the method followed can be 
applied to a more general class of surfaces of revolution. Further, it 
may be employed in the inverse problem of determining the form of 
the Earth’s surface, if the gravity be given. This problem was first 
solved by Stokest under certain hypotheses and by approximation. 
Helmert, Pizzetti, Almansi, and other authors have continued 
Stokes’s investigations. It appears that the method here employed 
is more general and better adapted to these investigations.§ 

1. I shall suppose an external Earth’s surface of equilibrium S to 
be a closed surface of revolution, having its axis of figure coinciding 
with the axis of rotation. Let the centre of S be taken for origin ; 

* See Stokes, ‘On attractions, and on Clairaut’s theorem’, Cambridge and 
Dublin Math. Journal, 4. 194, May and November 1849. 

+ See Mineo, ‘Sullo sviluppo rigoroso in serie di funzioni sferiche del poten- 
ziale esterno e della gravita superficiale d’un pianeta’, Boll. dell’ Unione 
Mat. Ital., Anno V (1926); and also Atti della Reale Accademia di Palermo, 
14 (1926-7). 

t Stokes, ‘On the variation of gravity at the surface of the Earth’, Trans. 
Cambridge Phil. Soc., 8. 672, April 1849. 

§ See Gulotta, ‘Sulla determinazione del Geoide, supposto di rotazione e 
di tipo ellissoidale, mediante le misure di gravita’, Atti della Reale Accademia 
di Palermo, 15 (1929). 
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let r be the radius vector, @ the angle between this radius and the 
axis of rotation. I shall take for the equation to the surface S 

r 2 | ' i 
_ 1+ (ay F yoy Pee? + (agg hot m2 ot mae et+. .. 
+ (Hn Fot+Ontot+- . -+Otng ls, e+. sing 


4 o i=j 
-=l4 5 oP ye) (1) 

j=1 i=0 
where P,; is the Legendre function of degree 2% in cos@; a, € are 
positive constants, « being less than 1; «;; are numerical coefficients. 
For particular values of «,;, equation (1) represents an oblate ellipsoid 

of revolution, a being the equatorial radius and « the ellipticity. 
Since Atk 
nek = 2 CF, 

i=0 
where Cp, ..., C,,, are constants, it then follows that every power of 
y admits of an expansion of the same type (1). We shall generally put 


-o § 


1 1 <j 

_— 2. (m)P__ <2 

a= onl! t 2 2% Pyje*!), (3) 
j=1 1=0 


for every positive integer m. 

Let M be the mass of the Earth and w the angular velocity. The 
potential V of the attraction for every point external to S is a 
spherical harmonic, which vanishes at an infinite distance, and at 
the surface S is equal to 

2 
F (6) =c—=r*sin®, 


c being a constant and / the constant of gravitation. 
Since sin?@ = 3(1—P,) 
sin*? = §(1—J5), 
we get, using (2), 


F(0)=c— 


1 


waz j 


af (1—P,+ > 
3f j= i= 
where the coefficients £;, are functions of «;;. 

It is easy to solve this particular Dirichlet’s problem. We shall 
prove that V can be expanded in the form 


Y 
=0 


BijPxe*), (4) 


; ae P, P.. . 
J rns ie sete: (5) 


r2 
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= AytAyet*+Aygett+ ..... phir 
Age? +Age*+ .. | 


. . . . . . 


ail 
A ae 24+ A n Fs + pd 


. . . . . . 


A,, being coefficients which can be successively evaluated. 

To prove this, we shall proceed by induction. Let V, denote the 
value of V at the surface S and suppose the equation 

V, = F(0) (7) 
to be satisfied as far as the terms in <2"-*; that is, we suppose that 
the }n(n+-1) constants 
Brg Aig, ++» Migs Bag, Bg 00es Beige 20s ee 
have already been calculated. We must then satisfy equation (7) as 
far as the terms in ¢?”; that is to say, we must calculate the n+1 
constants 
A,,, A», pee Ann» Ansin 

C = Cot C,e?+Cye4+...+¢,€%"+.... (8) 


Equation (7) may be written, to the desired approximation, as 


Let us put 


M n J : aes RE : 
— (14+ > > ofPPreY) + {Ai(1+ > > afPPaic) + 
j=1 i=0 


a j=1 i=0 


n—-1 j P. 
! 2 3)Pp _2j 2 2 
+Ane(1+ F > alPPye*)-+...+Arne™} $+ 


j=1 i=0 


| 2n Pyne 
+rAy es, n€ q2nts 


. wa? n j+i 
= Cy teyet +... p eye — A {1—P,+ > > By Pre? (9) 
3f j=1 i=0 

Since equation (9), by virtue of the preceding approximation, is 
satisfied as far as the terms in «?”-*, we must now equate the coeffi- 
cient of «?” in the two members of (9). But this equation contains 
Py, Pa, ..- , Ponsa only, and we have therefore as many equations as 
there are unknowns A,,, ..., 4,41,» C,- On equating the coefficients 
of P,, <2" in the two members of (9), we get an equation which con- 
tains only A,,,,, its coefficient being evidently 1/a***1. Our unknowns 

are then given at once by the equations 
1 


] 
aan = Gy, a4, — eer 


1 


q2nts3 


(10) 


n+l, n = On+1> 
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where o; (i= 1, 2,...,+1) are known quantities, linearly dependent 
on the coefficients last determined. 

By equating in (9) the coefficients of Pye*", we get 
Cn = %» 
which determines the constant ¢,,. 
It is therefore proved that equation (7) is formally satisfied by the 
series above proposed. 
2. Let v be the angle between the normal at any point of S, drawn 
outwards, and the radius vector; then 


cos v = Fs. dr\*\* sin an! enn 
"= Tab) J’  — 


On the other hand, we get from (1) 


a-°D > ir nae 


and since we have, as is well known, 


aP,, 
dé 
we find that w j 
cosv=1+e4 > > A,;Pe-* 
j=2 i=0 (11) 
~2 


= —{3P,+7P,+...+(4¢—1)P,-:}sin 8, 
3 ( 2i-1J 


oe 
sin vy = —e*sin@ ~ 2 Biz Po¢ re” 


Similarly, 1/ecosv admits an expansion of the same type as cosv. 


3. Let g be the force of gravity at any point of the surface S, then 
1 y 


cosy ér’ 
where W =fV +4wr* sin?6. 


We easily obtain the expansion of = by differentiation of (5), and, 
ér 


by using the formulae (11), we finally get for g an expansion of the 
form 


_fM 


a2 


-; wta+(Fu wa+3- ant 3 : vig Poie”. (12) 


4, Equation (12) gives g as a series of cecil s functions of 
even degrees and consequently as a series proceeding according to 
ascending powers of cos*#. It will now be shown that cos*@ can be 
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expanded in powers of sin®¢, ¢ being the latitude (that is, the com- 
plement of the angle between the normal at the surface S and the 
axis of rotation). 
First we have os 4 = cosvsin 6 — sin vcos 0. 
By means of (11) and (2), taking account of the relation P, = cos 8, 
we obtain 
vay~ —- f ¢ 21 Ee ee ae 
cos ¢ = sin O{1 + juy,(1+- 2P,)e? + [ (Apa t+ e412 — 3820) Ft 
2 3 4 41 
+ (Aya + $12 t+ Fo22) P+ (Azo t+ Fee2) Py let +. - .}. 
Squaring and remembering that sin?6 = 3(1—P,), it is easy to deduce 
an — for sin*¢ of the type 
sin?d = 4(1+-2P,)+ (09, Po +011 Po+ oo, Py)e?+ 
+(O92Po+042P2+ F22hy+ ogePo)et +. . + 
| a) > | I > 2 | 
+ (Fon ot+ointst+- --TOn+1, nt on r2)€ Meee (13) 
where the right-hand member contains Legendre’s functions of even 
degree alone, and consequently is a series of powers of cos*@. On 
reversing this series, we obtain, in terms of the coefficients «;;, 
cos*@ = {1-+-6a,,€?-+ (6a4.— 15a9— 6x9, X41 +3007, ) e4+...}sin 26+ 
+{— 6a, ,€?-+ (6091043 + 50099 — 6ay9— 9302, )et+.. jsin4’g+ 
+-{(63a7, —35ag9)e*+-.. .}sin®f+.... 
By means of this, the right-hand member of (12) becomes a series 
of powers of sind. 
It may be convenient to eliminate M by introducing the mean 
density k of the Earth. If + denotes the volume enclosed by the 


surface S, we put Wf Be 


(14) 


and 7 is given by the relation 


= ga a {fr >sin 6 dé. 
0 


Remembering the properties of Legendre’s functions of integral 
order, we easily get 
T = dra*{1+ 3a, €?-+ 3(a2, + apot+}az, )et+..-}, 
and M is then obtained in terms of & and of the parameters of the 
surface S. 
If S is an oblate ellipsoid, we have 


a ee is 
aH“? ey we Se 120° 
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If we omit terms in e*, equation (5) becomes 
y — M {wra® (> <r) ie ng “tt 

r' | 3 3 14f 1176f ~ |r 
+{- Sr e+( ee ears 
Similarly, equation (14) becomes 
cos?@ = (1—2e?)sin®f + (2e?— 5e*)sin4d + 4e4sin®¢. 
Let g, be the gravity at the equator and write 

wa 
Je 
We obtain for the gravity g at the surface, on neglecting also terms 


m= 


in me‘, 


5 1 3. ae Se ee. . 
~m—-—2— met) sin*¢ + » me? — — 4 sin‘. 
2 7 4 8 








INFINITE INTEGRALS IN THE THEORY OF 
BESSEL FUNCTIONS 
By A. L. DIXON and W. L. FERRAR 
[Received 6 May 1930] 


1.1. THE primary object of the paper is to prove many known, and 
some new, theorems about Bessel functions in a manner which is at 
once elementary and direct. Cauchy’s theorem on contour integra- 
tion leads, as we show, to a straightforward evaluation of such 
integrals as 


2) oO 


| J,,(ax)J, (bx) nay | J,,(ax)Y,(b2) ie 


x xr 


0 0 
and brings out clearly the discontinuity in the formulae expressing 
their values. 

Further, our method enables us to evaluate similar integrals in 
which other Bessel functions,* for example —Y,(x)+(2/7)K,(2), 
take the place of J,(x) or Y,(x), and to obtain a new, and in some 
measure elementary, proof of Nicholson’s integral formula for 
Tix) + Y7(x). 

1.2. The paper is a commentary on some parts of that wonderful 
chapter on infinite integrals, Chapter XIII of Watson’s Theory of 
Bessel Functions,+ and we use Watson’s notation throughout our 
work. In particular, we follow his practice (derived from Lommel) 
of using the symbols p, v to denote unrestricted numbers and reserv- 
ing m, n to denote integers. 

We first apply our method to a special case, where we write out 
formulae in full and develop the analysis in some detail. In later 
sections, we assume a wider knowledge of the properties of Bessel 
functions and omit certain details of the analysis. 


* Special cases of integrals involving L(x) = —Y, (x)—(2 7 )K (2) and 
M (x) = —Y,(x)+(2/7)K, (x) have been considered by Hardy, Messenger of 
Mathematics, 57 (1927), 113-20. 

+ G.N. Watson, A Treatise on the Theory of Bessel Functions (Camb. Univ. 
Press, 1922). References to this work will be cited as W. followed by the 
appropriate page reference. 
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Elementary example of the method 


2.1. Preliminary expansions. The expansions, near the origin, of 
Bessel functions of orders 0 and 1 may be written, for |z| = 3, 


Jolz) =1+0 (8), (2) = ie +0 (8°), 


2 
Yo(2) = = [log(Jz)-+-y]+-0 (8*log8), ¥;(2) = —— +0 (logs), 


7 
i = 1 
K 6(z) = —[log(4z)+7]+0 (810g), K,(z) => +0 (Blogs). 


The asymptotic expansions for large values of |z|—= R are con- 
tained in the general formulae, all valid in ranges which include 
—7 <argz<7, 
2\t 
H'(2) = Se) +i¥,(2)~ (2) ef" 1°01 4.0 (RY), 
mz 


H®(2) = Sse) ite) (Z) ete tT 14-0(R], 
TT: 


K,(z)~ ($) +0). 


When we try to discuss integrals like 
Jo(2)¥,(z) dz, f Jo(z)¥o(z) dz, 
0 
we see that the first is divergent at the origin. Accordingly, we 
introduce the function y,(z) to represent Y,(z) deprived of those 
terms which would give rise to divergent integrals if the series for 
Y ,(z) were integrated between the limits 0 and 1: thus, 


vile) =Y@)+—. (1) 


The functions h(z), k,(z), ete., are related in the same way to the 
functions H(z), K,(z), ete. 

2.2. ‘Increasing the argument by }3n.’ 

Consider the integral of H)(az)h'P(bz), where a and 6 are positive, 
taken around a contour composed of 

(i) y, the positive quadrant of a small circle | 

(ii) I, the positive quadrant of a large circle |z| 

(iii) the parts of the axes cut off by y and I’. 
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Using the formulae of 2.1, we see that 


HW (az)hP(bz) dz = O (8 log 8)? = 0 (1), 


Y 
$7 


47 
. . (at+b)Rsin 8 y —te—aRsin 
f H®(az)h(bz) dz <M J eatoyRsind dg +N [ R-te-aksind 9, 
r 0 0 
where M, N, for R > R,, are finite positive constants independent 
of R. Accordingly, since the integrand has no poles, 


R R 
i} HY) (ax)hP (bx) dx = + H)(axt)h (bat) dx +-0 (1), 
C) 65 


or, since the integrals converge when 60, Roo, 
f HW (ax)hD(bx) dx =i J H(axi)h)(bxi) de. 
0 0 


The process resulting in an equation like (2) is fundamental in 
most of the following work. The phrase ‘increasing the argument by 
17’ will be used to denote the complete process and, in general, the 
details of the process will be omitted. In some examples the equation 
corresponding to (2) will contain terms arising from the non-zero 
limits of the integrals taken over y and [. 

2.3. Formulae obtained by increasing the argument by 3x. 

We now transform (2) by using the formula 

Lori) (21) = e-”**K (2), 
together with its obvious corollary 
darth) (21) = e- 7", (2), 


and obtain, when a> 0 and b> 0, 
[ AP Cex Pox) dz = -sf K,(ax)k,(ba) da. (3) 
0 0 


A second formula obtained in the same way is 


2) 2) io 8) 


[ Ko(ax)h (bar) dx = i | Kg(axi)hP(bai) da = — | HS)(ax)k, (bx) dex. (4) 
0 0 0 

2.31. First case. a>b>0. 

For points z on I, since J,(z) = }{H)(z)+ H9(z)}, 


|H\D(az)J,(bz)| << M R-le-(a-b Resin 9, 
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where M, for R > Ro, is independent of R and 6. Hence, provided 
that a > 6, the integral along of H‘))(az)J,(bz) tends to zero as R 
tends to infinity, and so increasing the argument by 47 gives 


f H(ax)J, (bx) da =i i} H(axi)J,(bxi) dx = =| K,(ax)I,(bx) dx. (5) 
7 
0 0 0 


2.32. Second case. O0<a<b. 
Since we can no longer use the integrand of 2.31, we use 


J,(ax)H (bx) 


to make the integral over ' tend to zero.* Increasing the argument 
by 37, we get 
wn , 9; 
| )Jo(ax) HQ (ba) + ocd I,(ax)K, (ba) dz = ; +0 (8 log), 
: 7 b 
ra) 


where the terms on the right represent the integral over y. Hence, 


q én 
| Jg(aee) HY (bx) += To(azr)K, (bee) dar = ; (6) 
{ 

2.33. Third case. a=b>0. 
On increasing the argument by $7 with H )(ax)J,(ax) as the 
integrand, the integral along I will tend to a non-zero limit, namely, 

ia 
: 9 1 - 1 
am. f SH (az)HQ (az) dz = | ili d) = — a’ 
r 0 
and it is this term, zero when the integrand is H)(ax)J,(bx) and 
a >b, which marks the discontinuity of the integral of J(ax)J,(bz). 
Hence (5) is replaced by 


{ 


0 as 
)} H\))(ax)J,(ax)— = K,(ax)1,(ax)} dx = = (7) 


while (6) is replaced by 
rt 25 
f | Jo(aae) HP (axe) + = I,(ax)K,(ax)\ de = 


0 


1 
—. (7a) 
2a 
* Note that h(bx), while avoiding difficulty at the origin, would make the 
integrand on T tend to infinity as R increased. 
3695 K 
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2.4. Hardy’s theorem on the integral of M,(ax)L,(bx). 
Taking real and imaginary parts in (3) and (4), we have: 


Whena>0,b>0. 


wo 


fi J,(ax)J,(bx)—Y,(ax)y,(bx)} dx = — S| Kolax (02) dx, 
7 


fi {Jo(ac)y, (ba) +-Y,(ax)J,(ba)} dx = 0 
f Kat (ax) J, (bx) dx = — fn ax)k,(bx) dx, 


[ Ko(ae)y (bx) dz = f Y,(ax)k, (bx) da. 
0 0 
But, as the expansions of Y,(x) and K,(x) show, 
9 9 
Y;(x)+ — k(x) = Y,(x)+ - Ky (2), 
7 7 
so that we can combine formulae (8)—(11) in such a way that y,(x) 
and k,(x) are replaced by Y,(x) and K,(x). Using the notation 
adopted by Hardy,* namely, 


L,(«)= —¥,e)—=K(z), Ma) = —¥,(@)+ 


v 


we obtain from (8) and (11) 


fi J,(ax)J,(bx)+-Yy(ax)L,(bx)} dx = = | Ka(ax) La(b2) dx, 
7 
0 


and so, since the integral of J,(az)J,(bx) is convergent, 


| Jo(ax).J (be) dx = f M ,(ax)L, (bx) dex. 
0 0 


From (9) and (10) we have 


| (lb) Lo(ax)—Jg(asr)y, (b)} dz = =| J,(ax)k,(ba) dx, 
0 
which may be written as 
| { Jo(ax) L, (bx) + L,(ax)J,(bx)} da = 0, (14) 
0 
where the integral of each term is — provided that a +b. 
* Messenger of Math., 57 (1927), 113-20. 
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2.41. Whena>b>0. 
Taking the real parts of (5) and (6) gives at once 
oo ao 1 
f Jo(ax)J, (bx) da = 0, f J(ax)Jo(bx) dx =—, 
0 0 
which, when combined with (13) above, gives Hardy’s results 
i) eo 1 
J Mo(ax)L, (be) dx = 0, f Zax) M (bx) dz = => 
0 0 
Taking the imaginary parts of (5) and (6) gives 


oo 
ine) 


J Yala) (ber) die = = | Kalan), (bx) dz, (17) 


r r 
f {¥, (ax) Jo(bee)-+ = Ky (asx) Io(be)} dx = 0. (18) 
0 
In (17) we may substitute the series for /,(bx) and integrate term 
by term using Heaviside’s formula 
f K,(t)t# dt = 2#-*AT(4y)}?. 
0 


We then have 


i) oO 9 ° 
‘a f Jo(ax)y, (ba) da = J Y,(ax)J, (bx) da = = | K (ax) 1, (bx) dx 
7 
0 0 0 
l ~ l b\2mt+1 
= — nd 19 
7a pS m+1 (5) (19) 
which, for comparison with less special results, may be written in 
the forms 


1 a* b . 
— log —, pe? waathi(1, 1; 255). 


ab a*— ss 


Again, on subtracting (18) with 6 = 0 from (18) with b 4 0, we have 


0 


J Y, (ax){Jo(bx:)—1} dx = — = | Ku(axy{ (be) 1} dx 
7 
0 


a*—b? 


1 
| pee 
—s 


K2 
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2.42. Whena=b>0. 
The real and imaginary parts of (7) and (7a) give the formulae 


f ae ax)L,(ax) dx = f foo \,(ax) dx =>, 
2a 


i! Y,(ax)J,(ax)— 2 K,(ax )Iy(ax)| dz = 0, 
a 7 


0 
a) 


[ {laa (ax)+ - . ~T, o(ae)K,(ax)} dx = 0. 


ty 


0 
The integrals in (22) and (23) are convergent, but the integrals of 
the separate terms are divergent. The dominant term at infinity 
of the integrand in (22) is, for example, 
(1—sin 2ax)— ; 
Tax TAX 


The discontinuous integral of Weber and Schafheitlin 
The integral associated with the names of Weber and Schaf- 


heitlin is 00 


F (a x) J, (ba) 
‘< dx. 


Our method of ovebaiiiin’ is relatively simple. Having put suitable 
restrictions upon the values of A, , v, we transform the integral by 
means of Cauchy’s theorem [see (28) and (29)] and find at once the 
discontinuity in form which marks the two cases a>b and a <b. 
We then expand the new integrand in a power series and integrate 
term by term, a procedure which is readily justified. Finally, analytic 
continuation, applied to the integral gua function of » and v, dispenses 
with unnecessary restrictions upon the values of these variables. 

3.1. The standard integrals. 

All the integrals which occur in this section depend, for their 
evaluation, os the two formulae* 
AHP (yt yf pax i (bn) dx 

a 
0 


ba bck b 
= q\--1h",F, pent a = 4 I By ne 


s+ 
> >? 


2 ~- 


xT (Se \n(- re) 
2 2 


* Neither result is new; (24) is an obvious modification of W. p. 410, 13.45; 
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K,,(ax)K (bz) i 
aa 


aera—2)| 


0 


eee 2 a? 


- 


xr (bette) p(t ae), 


ee —_— 2 
= a-o-tbreF(? Atuty 1—A err. 1—); 1-3) x 


_ — - 


The first of these is proved by substituting for J,(bx) its series 
expansion and using the result* 


J K (8 dt = 28 (F)r(* 5") 


to effect the term-by-term integration. 


The second is the result of writing 


K fia) = in 1_,(ba)— T (bx) 

sin vz 

and combining the two hypergeometric series, which arise on using 
(24), into a single hypergeometric series. The formula used to effect 
the combination is given in Whittaker and Watson, Modern Analysis 
(1920), p. 291, and, to apply this formula in the present instance, 


all we need do is to write 
h=l-y, pa—P, v=atB—y, 2=(a*—b*/a%, 


In both formulae we suppose the values of a, b, A, w, v to be such 
that the integrals converge. In the second formula we supposet 
either that 0 <b=a or that a> 0 while |b] <a and |argb| < $n; if 
neither condition is satisfied, the relation we have used between the 
three hypergeometric series is no longer valid. 


(25) is given by Titchmarsh, Journal London Math. Soc., 2 (1927), 98, who 
proves it in another way for real values of a and b. 
* W. p. 388 (8). This formula is proved at once by substituting for K,(t) 
; _ 14). T . 
its equivalent integral representation (3) -G¢) | e~(u® — 1)”-4 du, and then 
I'(v+ 4) 


integrating first with respect tot. The corresponding formulae for i) J;,(t) tel dt 
0 


ey 
0 


and | ¥,(t) ##~'dt are derived from it by ‘increasing the argument by 37’; 


they are particular cases, when b = 0, of (28) and (30). 
+ If a is not real and positive we can easily transform and make it so. 
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3.2. Evaluation of the Weber-Schafheitlin integral. 
Suppose that a >b > 0, and, to ensure convergence in what fol- 
lows, that A, 1, v are subject to the conditions 


RQ)>—-1, R(+p+v+1—d)>0. 
When we take as integrand H‘)(az) J,(bz)z> and increase the argu- 
ment by $7, we obtain 


[ Helen) ae HY(axi)J, (bei) 
x xa 

0 0 

*K (ax)I,(bx) 


ar 


9 
sail = eh(v—p—Ari Led 


7 


dx. 

0 
Taking now as integrand H®(az) J,(bz)z> and this time decreasing* 
the argument by 37, we obtain 


ee ree =e-Hr-n-bet | K ,(ax)1,(ba) nF (27) 
xn 7 xn 
0 0 


Adding (26) and (27), we have 
oe) 2) 
dx = —cos 


x 7 “ 
0 0 


a ———— 9 lk A. — dz, (28) 


xr 


ora) 4) 


te »; {2 (ax) J,(bx) 2 a —_ dx, (29) 


dx = —cos 
x = x 


0 0 


the second form being an obvious modification of the first. 

We have thus evaluated the Weber-Schafheitlin integral when 
a+b. The discontinuity in the formulae is apparent. In order to 
write the formula in the shape given in W. p. 401, we must express 
cos $(A+-—v)z in terms of gamma functions and use (24): the result 
is, when a >b, 
ba T(—+it+hth) op 
DY al A AY Sof) 

AD (y+ IP ($A+ du— v4 4) 
and is valid provided that the integral on the left of (28) is con- 
vergent, i.e. that R(A) > —1 and R(u+v—A+1) > 0. 

* That is, the contour used is the reflection in the real axis of the previous 

contour. We now use HY, (axe~*7") = - eh! (axe*™"), cf. W. p. 75 (6). 
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On subtracting (27) from (26) we have, when a > b,* 


x 
xr 7 2 
0 0 


— ee Bin OMe | Kyla) dx. (30) 
x 


3.3. The critical case a= b. 
If = R(A)>0 and R(u+v—A+1) > 0, then for }<a<1, 


[7 OS, (at) | O(x-, 
tr T 
x 
and it is easy to show that 


: 2 


"J, (at). . 
(at) J, (at) on | — di 


tA b= 
0 


- d-TAT(—P+e+b+4) | 
QT(3A+ du—dv+ 4) (GA4 3u4 4-2) (GA— 3-49) 


Similarly, the critical case of (30) gives 


(31) 


fA 


Y,,(at).J,(at) cies —_— di (32) 


b—a 
0 


> 


sin }(v—p—A)a a PAjr(—3 sA— -de+4 5+ 3) I'(— A+ Sut 4v+4) 
. 7 an Sut 4v+3)0(SA—3u4 +3) 
provided that R(A) > 0 and R(+y+v—A+1)>0. 
3.31. When p—v ts an odd integer. 
Changing the notation, we consider 


[ i —p- _;(at) dt, (33) 


0 
where p is a positive integer or zero. The asymptotic values of the 
integrand show that the integral converges at the upper limit not 
only for R(A)>0, but also for R(A)>—1. There is no need to 
calculate the values afresh as these can be derived from the values 
obtained under the more restricted conditions of 3.3. 


* The formulae for J J,,(ax) ¥, (bax dx and f Y,,(ax) ¥, (bx)a4 dx appear 
later in (39) and (40). 
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For any « whose real part is positive, (31) gives,* when 
0 < R(A) < R(2a), 


i 8) 


Tu. p(Mt) Jp (at) di 


tr APL p+ ore 


| 


IP(A)P'(a— $A) 
LA) P(a+)P(A—p)’ 
0 
Now the integral converges at the lower limit if R(2«—A) > 0, and 
may be replaced by 
(+0) 
1 | Jy (at) Ip (at) di 
e27i(2a—-A-1)__] tA ‘ 

provided that 2~a—A #n, a positive integer. In this form it is easy 
to prove that, for a fixed «, the integral (33) is an analytic function 
of A in the strip R(2«) > R(A) > —1 when points A= 2a—n are 
excluded, and is continuous at the excluded points. Hence if R(a) > 0, 
(31a) is true for R(A) > —1 and, in particular, 
i T(1—}A+p)sin(jA—p)r —_ (—)P 
Jv p(t) Jy» (at) dt = lim Sr aes ( = ; 
(31b) 

Similarly, by taking a fixed A which satisfies the condition 
R(A) > —1, we can show that (31a) is valid when R(2«) > R(A) > —1. 

On comparing (31a) and (31 b) with the values obtained by putting 
w=atp, v= a—p—l1 in (28), we see that 


(i) when R(A)> 0, the expression for foto (at) Jy _»—y(bt)t dt 


changes its shape as the ratio (b/a) increases through unity, but, 
considered as a function of (b/a), the integral is continuous, 

(ii) when 0 > R(A) > —1, the formula changes and, moreover, the 
integral is discontinuous. 

3.32. When p—v is an even integer. 

By the method of 3.31 we obtain from (32) the resultst 


J, ,»(at)Y,_,,(at) Pe a\—sin(p—}A)n T(p+4h— 
tr ~ 2QAD(1—A)sin Ax © (p+ 


0 


* W. p. 404 (3) has a curious misprint, }A+-} for A, in giving this formula. 

+ The ge ¥ = an integer must be excluded at first. But since 
Y,(z) — Yn(z) uniformly with regard to z, except near z = 0, the integral 
values of « are easily included in the final results. 
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valid when p is an integer greater than }R(A—1) and 
R(2a+1) > RA) > —1, 


5 


(— 
and | J+ (M)¥ ~p(at) dt a ~ a 


0 
valid when p is zero or a positive integer and R(«) > —}. 
3.4. Some associated formulae. 
If a and b are positive and we assume, in order to ensure con- 
vergence, that A, », v satisfy the conditions* 
a=b; 0<RA)<1; — |Riy)|, |R@)| <3—}RQ), 
aA#b; —1<RA)<1; = [Riy)|, |Rr)| <}—}/RQA)|, 


increase of the argument by 47 gives, whatever constant values 
6 and ¢ may have, 


Ky, (ax)+ 0K ,(—aaxt)}{K, (bx) +K,(—bzxi)}a- dx 


- ie [ {K., (avi) + OK ,(ax)}{K, (bat) +-$K, (bx )ya-A dav. 
0 


This gives, on writing 
K (axi) = —4aie-¥ "HH (az), 
K,(—aaxi) = 4rie?™ H™ (az), 
the following formulae: 


ie 2) 


(ax) i (bx) die nit Sm (A dx, 
e H 


0 
; fH) (ax) (ber) a “oem (Sore 
an 


- — é2, 


nH 
{ 


x 
) 0 


a i) 


ax)H (ba) ae wie Wenvont (i (ax yee (bx) | dx (37) 
; x 


0 


eK (ax\H (ba ¢ Wa) r 
K,,(aa aa (bx) dx = jlsnvomt | a = y(bx) i (38) 
x x 


0 0 


* In a later section, 3.6, we show how our results can be applied when some 
of these restrictions, those needed for convergence at the origin, are removed. 
+ The dominant term in the asymptotic expansion of Ky(z) is (7/2z)te—. 
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From (35) and (36) we have 
y (ax) Y, (ba) 


Be 
uA 


dx 


= dx, 


( ¥,,(ax)J_,(bx)-+ (4/m2)cos }(A+-p-+-r)K,,,(ax)K (be) 


. dx, (40) 
x 


0 
and we see, on referring back to (28) and (30), that we have now 
evaluated the integrals which arise from all possible combinations of 
J and Y.* 
Q5 7 Y ty, > . . 
3.5. The functions L(x), M(x). 
If in (34) we put the particular values 
—byni — —bvri 
6 = —e-te7, b = —e-t7i, 


we obtaint 


fi L,,(ax)-+-tJ,,(ax)}{ L, (bx) +-tJ,(ba)}a- A da 
= e~hA+h iat (D(a) — iJ, (ax) }{ L, (bx) — iJ,(bx)}a* da. 


Write now «a= }(A+p+v—1)z. Then (42) shows that 


ext [ L,,(ax)+iJ,(ax)}{ L,(bx)+%J,(ba)}a~ dx 


is unaltered when 7 is replaced by —?. Accordingly, 
8) 
cos x {LZ,, (ase), (bir) + L,(ba)¥J,,(aa)}a- dx 


0 
a) 


— —sin a i} {L, (ax) L,(bx)—J,,(ax)J,(ba)}a- da. (48) 


0 


* It is, of course, possible to evaluate (30), (39), and (40) by expressing Y 
in terms of J, but (39) in particular becomes tedious when thus attempted. 
{ The functions Ly(x), My(x) were defined in (12). 
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From this it follows that 


x 


res (ax)sin « + J,,(ax)cos a}{ L,(bx)sin « + J,(bx)cos aja-A dx 


pL 
= (cos*a + sin*a) { J,(ax)J,(ba)e dx. (44) 
0 


Changing the signs of both @ and ¢ in (41), we prove similarly that 
f {M,,(ax)sin « + J,,(ax)cos a}{ M,(bx)sin « + J,(bx)cos ala da 


= fJ,(ax)J,(bx)x> de. (45) 


Again, changing the sign of ¢ only in (41), we see that 
ef L,(az) +iJ,,(ax)}{ M,(bx)+iJ,(bx)}a dx 
changes sign em i is replaced by —1. Accordingly, 
cos x { Lyla) M, (bx) —J, cer), (b2)}e- dx 
0 


= sin a J {L,(ax)J, (be) + J,,(ax)M,(ba)}e dx, (46) 
> a 
and from this it follows that 


| (L,,(ax)cos a — J,,(ax)sin «}{M,(ba)cos « — J,(ba)sin aja da 


= (sin’a + costa) { J, (aix)J, (bare dx. (47) 
0 


Moreover, with the aid of formulae already proved, (43)—(47) pro- 
vide simple machinery for evaluating, in terms of hypergeometric 
functions, such integrals as 


J E.lax)L, (be) dx, [ L,ax)M, (bx) dx. 
0 0 


3.6. The use of circuit integrals and other modifications. 

All the formulae of this section have been obtained under restric- 
tions on the values of » and v which ensure that the integrals con- 
sidered shall be convergent at the origin. The corresponding formulae, 
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when these restrictions are removed, can, in general, be obtained by 
substituting for the line integrals from 0 to 0 a proper multiple of 
the circuit integral which starts from 00, makes a circuit round the 
origin, and returns to 00.* The evaluation ‘n terms of hypergeometric 
functions will, in general, be unaltered in form. 
Also, this circuit integral can always be replaced by a line integral 
with a suitably chosen integrand.+ For, consider the integral 
(0+) 
f x?-lf(x) dx, 
% 
in which f(x) is regular at the origin, and p lies between the negative 
integers —n and —n—1. It is equal to 


(0+) 


| UI fle) —f(0)—af (0)... — Ef (0)} da, 


that is, to 


© 
(e27tp — 1) [ar } f(x)—f(0)—af'(0)—. — =f (0)} dx ; 
J n! 
so that, except when p is an integer, we can substitute this line 
integral for the circuit integral. 

Another method by which restrictions on the values of » and v can 
be avoided, namely, the use of y,,, k,, for Y,, K,,, has been exemplified 
in § 2, where we considered the simple case » = 1, v= 0. The results 
of 3.5 in particular, when modified by this method, yield some 
interesting formulae. 

3.61. Formulae when p and v are integers. 

Following the elementary case of § 2, we write 

<— (—)"(2n—m—1)! 


m!(4z)?" 2m , 


m 0 


n 


z —)"(2n—m)! 
= K,,,,(z)—}4 z i‘ 
) 2n+1( ) 2 : m!(4z)2"+1 2m ’ 


m 

and use 1,,(z), m,(z) to denote similar modifications of L,,(z), M,,(z). 

It then follows directly from 3.5, when we consider the two cases 
sina = 0, cosa = 0, respectively, that 
* Cf. the use of such an integral in 3.31. 


7 Cf. A. L. Dixon, Messenger of Math. 33 (1904), 176. The method is that 
used by Cauchy, and later by Saalschiitz, in the theory of the Gamma function. 
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(i) when r and s are positive integers and r-+-s is of the form 4N-+1, 

J L(ax)m,(ba) da = J m,(ax)l,(ba) dx = f J,(ax)J,(bx) dx, (48) 
0 0 


0 


x 


{l.(ax)J,(ba)+-J,(ax)l,(bx)} da 





= J fom, (ax)J, (bx) +-J.(ax)m,(bx)} dx =0. (49) 


(ii) when r and 8 are positive integers and r-+-s is of the form 4N-+3, 


io 2) 


[t(azyi,(bx) dz = [ m,(acx)m,(ber) dz = f J:(ax)J,(b2) dx, (50) 


f E(ax)J, (bar) +-J,(aar)m,(bar)} dx = 0. (51) 
0 
Finally, it may be pointed out that 1,(x)= L,(x) when r = 0, 1, 2 and 
that m,(x)= M,(x), so that the above include Hardy’s formulae* 
for f L,(ax)M,(ba) dx and f L,(ax)L,(bx) dx. 


Nicholson’s Integral} 


In this section we prove Nicholson’s integral formula, 
eo 
J?(a)+ Y¥2(x) = = | K,(2x sinh t)cosh 2vt dt, 
0 
by means of more elementary Bessel function formulae. One of these, 
namely (52) below, is known to be valid when |argz| < }7; since we 
use it when argz = +47, we begin with a preliminary theorem which 
justifies such a procedure. The proof of Nicholson’s formula starts 
with 4.2. 
4.1. A preliminary theorem. 
On p. 181, 6.22 (8), Watson gives the formula 


K,(z) = | e-*°*"! cosh vt dt, 


og 


* Hardy, loc. cit., evaluates the integrals by another method. 
+ W. p. 441, 13.73. 
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valid when |argz| < 47; and on p. 180, 6.21 (11) and (10), 

K (ix) = —4hrie-#™'H®)(x) = | e-** cosh cosh vt dt, 

0 
K,(—iax) = drie™ H(z) = feist cosh vt dt, 
0 
valid when x > 0 and —1 < R(v) < 1; so that, when —1 < R(v) <1, 
the first formula is valid for jargz| < }z. 
It will follow that the formulae* 


. ft) 2[ Kysy(2z cosh t)cosh(u—v)t dt 
0 


“a 2( K,-,(2 cosh t)cosh(u-+-v)t dt 
0 


are valid for |a1 
when |R(u-+v)!| 


‘'gz| = 4n, the first when |R(u—v)| <3, the second 
<j. The proof follows the lines of that given by 
Watson, but, since the integrals involved are not absolutely con- 
vergent, additional care is needed in handling the transforma- 
tions. 
4.11. Proof of (52). Suppose first that | R(u)| <4 and |R(v)| < 4. 
When z is real the above formulae give 
K,, (ta) K, (iz) 
R R 


= lim } f e708! cosh pt dt x e~tz cosh“ cosh vu du 
R->@ 
i —-R 


= lim i[jenem +eosh“) cosh yt cosh vu dudt, 
R>& 

Sr 
where the double integral is taken over the square Sp bounded by 
u=+R,t=+R. ; 

We can change the variables from ¢, u to 7’, U where t+-u= 2T, 
t—u = 2U and consider the integral to be taken over the square Qp 
bounded by U = +R, T = + BR if we can show that the integral over 
such a triangle as that bounded by t= 0, u= R, t+-u = 2R tends to 


* W. p. 440, 13.72, where the formulae are proved with |arg z| < }7 for 
any p and vp. 
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zero as R tends to infinity. The hypotheses | R(u)| <4, |R(v)| <4 
ensure that this integral does tend to zero ; for then, 


e- iz coshu cosh vu du (53) 


2R-t 


22 2° 
aii =| «teem , 2 cosh vu ie. | en ix cosh u d 7( sinks) du 


xsinhu |p x du\ sinh u 


= 0 (e*®)+0 (e+), 
and the double integral over the triangle is 
R 
fe (e444) dt =o(1). 
0 
We have then, when z is real, 
K (ix) K, (tx) 
RR 
= 2lim f fe cosh Teosh U cosh(u+-v)T7' cosh(u—v)U dTdU. 
R> 
00 
If now we restrict » and » still further and suppose that | R()| < }, 
|R(v)| <4, then | R(u+yv)| <4, and we can prove in the manner of 


(53) that 


fe-sircom T cosh U cosh(u+-v)T dT =o (e-?*), 


R 
R 
ff e-*iz cosh T cosh U eosh(u+-v)T' cosh(u—v)U dTdU 


0 


R 
= [0(e#)|cosh(u—v)U | dU =o0(1), 
0 


Accordingly, we may write 
K (ta) K, (tx) 


pe 
Raw 


= 2 lim | eeaaas Tcosh U eosh(u+-v)T cosh(u—v)U dTdU 


R—> 
wae 


= 2 [ Kys (dia coshU )cosh(u—v)U dU. OG 
0 


But, by using the procedure of (53) again, we see that the last integral 
converges for real values of x whenever | R(u—v)| <3, and the theory 
of analytic continuation enables us to assert that (54) is true pro- 
vided that this condition holds. Similarly, the second form of (52) 
is proved to be valid when | R(u+-v)| <3. 
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4.2. Increasing the argument by 4x. 

It is now convenient to start from the other end and write, with 
R(x) > 0 and |R(a)| <1, 

[Ka(2e sinh t)e” dt = | K ,(2au){v(u?+1 +up? 
0 

If further R(b) <3, the latter integral, when the argument is in- 
creased by $7, becomes*: 


du 
(u?-+1) 


: lu 
K,,(2txu){v(1—u?)+ iu}? : 
[Ra } v(1—u?) © 


du 


+ | K,(2ixu)fiv(u?—1)+iu}? ; 
J al Mt ( ) 1 } \(u?—1) 


In the first integral put uw = cosé@ and in the second u = cosh?t, and 
we obtain 
a2 
e-tni [ K,(2x sinh the dt 


0 
ba 


oe 
= if K ,(2ix cos0)e-% dé + [ Kal2ia cosh t)e” dt. (55) 
0 
Now put a= p+v, b= +(u—yv) and use (55) to transform (54). This 
gives, when 2 is real and positive, 


K,, (ta) K, (tx) = [ Kyusv(2ia cosh t){e"—"¥ + e-#-"™} dt 
0 


= | K,,,,(2asinh t){e*e—77eu—-l eter) rig—u-)} dt — 


pty 
— 2% fx, K,,.,(2ia cos@)cos(u—v)0 d0. (56) 


Multiply each side by e#+)7? “a add to the conjugate result 
obtained by writing —7i for 7; we thus obtain 


—4n*(J,(0)J, (x) -¥, (0)¥,(@)} 


an 
= 2/K,, 2x sinh t){e-”" cos vm + e+ cos rr} dt — 
0 


20 fy u+v(2% COS8)cos(u—v)O dé ; 


* The contour is to be suitably indented at the point u = 7. 
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or, using the fact* that the last term is equal to —7*J,,(x)J,(2), 
J,,(x)I,(x)+Y,()¥, (2) 


= 5| K,,;,(2x sinh t){e#—”¥ cos vz + e“#-”* cos pr} dt, (57) 


0 
the formula being valid} when R(x) > 0 and |R(u+v)| <1 
As we shall see in the next section, Nicholson’s formula is an 
algebraic variant of (57). 
21. Formulae derived from (57). 
Changing the signs of » and v in (57) gives 


K , (2a sinh t){e“"—-”™ cos va + e#-”* cos pz} dt 


J (2 v)J_ - y(a ss p(x)Y_ p(x x) 
(J, (2), (ae) +-¥,,(a)¥, (2e)}08(4—v)er + 

+{J,(2)¥,(x)—Y,, (x)J,(x)}sin(u—v)z. 
Hence, combining a with : 


Jule)¥(a)— Y,()J,(x) 


4 . : : . = 
——| K,.,(2asinht)fe-¢- sin ux — e#-”" sin vz} dt. 58 
= pty t B 

0 
Again, we have the two equations 


J_,Y,—Y_, J, = J. ¥,—Jd, Y, eos ua — (J, J, +Y, Y,)sin pr, 

J, Y. y—Y,J_,= Y, ¥,—d, ¥, eos vm + (J, J,+Y, ¥,)sin vz, 
sO that, on changing the sign of » alone, and then of v alone, in (58), 
we can derive the formulaet 


J,,(x)J,(%)+-¥, (x) ¥, (x) 


); 


io 3) 


5| K,_,(2x sinh t){e#+"¥+-e“#" cos(u—v)a} dt, (59) 


_4 mip | K,_,,(2x sinh t)e-“*”" dt, (60) 
m2 


0 
valid when |R(u—v)| <1 and R(x) > 0. 

* W. p. 150, 5.43 (1). 

+ The henctiaite in (57) are analytic in R(x) > 0 and have been proved 
identical along the positive real axis. 

t These are the formulae (4) and (5) of W. p 


. 445, 13.73. There is a misprint 
in Watson’s work, where e'“*”” and e ~uton 


should be interchanged. It has 
been noted by Hardy, Proc. London Math. Soc., 23 (1925), 1xi, who indicated 


independent proofs depending upon Mellin’s formula and Lerch’s theorem. 
3695 A 
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Putting » = v in (59), we obtain BBE 3 formula 


> r9 8 r . 
J?2(x)+ Y2(x) = 5| K, (2x sinh t)cosh 2vt dt, 
a 
0 
and, making »—>v in (60), we obtain Watson’s formula 


eY,(z)_y (x) ad.,(a) 


Ov Ov 


J,(x) 


v 


4 se ; 
—— | Kote sinh t)e-2” dt, 
7 


0 
both valid for all values of v if R(a) > 0. 
Making » —v in (57) and (58), we obtain 


8cosvz [ , ’ 
J?2(x)+ Y¥2(x) = a | K,,(2x sinh t) dt, 
pe 
0 
OY (x . J, (2 . ( 
J,(2) : v(a - Y (x) . (a i -4/ K,,(2x sinh t){7 cos v7 —2tsin v7} dt. 
Ov Ov _ 
0 
4.22. Formulae derived directly from (52). 
Writing z = iz and z = —iza, where « is real, in (52) we obtain, on 
adding and subtracting the results, 


J(u) J(u) —Y,,(@)¥,(x) = | Y,,.,(2a cosh t)cosh(u—v)t dt, (61) 


7 
.y(2a cosh t)cosh(u—v)t dt. (62) 


Putting » = v in these equations, 


| Y,,(2a cosh t) dt, (63) 
7 


0 
a) 


9 
—- | Jen(22 cosh ¢) dt, (64) 
7 
0 
the formulae being valid for all values of v. 
On changing the sign of » we get two formulae analogous to the 
results (61) and (62). Further, when we put » = v we obtain 


S(z) =x | tole" cosh t)cos vz —J,(2x cosh t)sin vz}cosh 2vt dt, 


‘ (65) 


= — 5] 0 Y,(2a cosh t)sin vz-+-J9(2a cosh t)cos yz}cosh 2vt dt. 
(66) 
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By applying the method used in (53) we see that these results hold 
whenever |R(v)| < ?. 
4.23. Complementary Formulae. 
The set of formulae, so far obtained, for evaluating the integrals 


( K ,(2(e-Le+))e dt can be completed, if we notice that, for a re- 
stricted range of values of a and b, pw and v, (55), (57), (58) are still 
true when we substitute + ix for x, that is, when argz = +432. 

The restrictions necessary to ensure the convergence of the in- 
tegrals involved, will be assumed made at each stage. 

Transformation of (55). 

We begin by substituting —ix for x in (55). We have 

eon | K ,(—2ta sinh i)e™ dt 


0 
) 


4 - 
--a | K (2a cos 0)e-% dé +] K (2x cosh t)e™ dt. 
0 0 
Changing the sign of 7 and subtracting, 
e-toni [ K,(—2ixsinh t)e” dt —e*' | K,,(2ixsinh t)e” de 
0 0 in 
= 23 | K,,(2x cos 6)cos b6 dé. (67) 
0 
3ut [W. p. 441, 13.72 (2)] 
$7 
ee Ln 
| K (2a cos 8)cos bé dé = ees {1 _ya-»! ya +b) ~~ Tya-w! aso} 


which is (in either form) an even function of b, so that we have also 
the identity, 


oa oe 
e~tont | K,(—2ixsinh t)e” dt —elbri | K,(2ie sinh t)e” dt 
r 0 


= elbri | K,,(—2ix sinh t)e- dt —e-¥7 [ K_ ,(2ix sinh t)e- dt, (68) 
0 0 
or, in the more usual notation, 
( J,,.+,(2x sinh t){e'*—™ cos vr —e”-#" cos ur} dt 


0 oo 


= | ¥,.,(2xsinht){e#-™ sin vz—e?-»"sin wx} dt. (69) 
0 


2 








144 A. L. DIXON AND W. L. FERRAR 


Transformation of (57)+-7(58). 
Combining pss and (58) we find that 


Hi@)HP@) = 24,2) HH ( 


4 , 
ce 4 9+ oj f ol Mo—vai_| olv—pihp pt 
= | Kyi,(2zsinht){er—e! + leu} dt, 
é 


0 


so that, on writing ia for z, we obtain 


9; 
4 ; , " 
—K pl) K, (x) + 2ee7'T (x) K, (x) 
i ) 
a | {J,,+y(2a sinh t)—tY,, ,,(2a sinh t)}{e"-™e-¥7! + lew} dé, (70) 
from which, after ee real and imaginary parts, 


2a sinh t){e’-/cos ya +- e'#—-cos vr} dt + 


2 cos pr ,,(: [4 u+v(2 


+ fY, Y,,,(2x sinh t){e?-»" sin ua — e#-” sin vz} dt, (71) 
0 


i 8) 
. 


5) 
— =cosprK ,(x)K,(x) = sin(u+v)z | J,,,,(2asinh the” dt + 


0 
+ -[Y, Y,,.,(2a sinh t){e“—” cos(u+-v)a+-e"-9 dt. 
0 
From (72), on interchanging » and v and adding, we get the sym- 
metrical formula 


= see k(u—v)n K ,,(«)K,(x) 


7 


= sin }(u+yv) f J,,.,(2x sinh t)cosh(—v)t dt + 


0 
an 


+cos $(u+v)z [ Y,,,(2” sinh t)cosh(u—v)t dt. (73) 
0 
Special cases of this are [uy =v; p = —yv] 
—7 sin vr | -J,,(2z sinh t) dt—7n cos vm [ Y,, (2% sinh t) dt, (74) 
0 0 


ao 
” 


= | Y,(2x sinh t)cosh 2v¢ dt. (75) 
COS V7T 
0 








THE THEORY OF BESSEL FUNCTIONS 
Also from (71), making use of the identity (69), we obtain 
eo 
= | J,.4,(2xsinh t)e”- dt, 
0 
valid if p+v > —1, v—p < 3; and, changing the sign of v, 


= | J,_,(2xsinh t)e-*! dt, (77) 
0 


valid if u—v > —1, —v—p <3. 


If in (76) we put »=v, we obtain the known formula [W. p. 435 
13.6 (3)] 
K, (ak), (Sak) = (eee (78) 
ue 1 


and from (77), we obtain 


K,(x) I(x) [ Jo(2x sinh t)e-* dt. 
0 











THE POTENTIAL OF AN ELECTRON IN A SPACE- 
TIME OF CONSTANT CURVATURE* 
By H. 8. RUSE 
[Received 14 May 1930] 


§1. Introduction 
THE present paper contains a solution of the equations of electro- 
magnetism as they existed before the introduction of Einstein’s 
Unified Field-Theory. 
Let X,, denote the electromagnetic six-vector in a space-time 
specified by 2 
ds? — Jyydxtda”, (u,v = 0, 1, 2, 3). 
Then it is well known that Maxwell’s equations may be reduced 
to the simple form 
Cb, od, 
pv a — ae ? (1.1) 
ox” = ath 
and Ae? = JP, (1.2) 


where ¢,, is the electromagnetic potential-vector, J’ is the charge- 


xX 


and-current vector, and X#” is the divergence of X¥”. 


Substituting from (1.1) in (1.2), we are quickly led to the familiar 


equations ~* ‘ 
9? by ep—($>)y+ Gide te J, (1.3) 


where G¢ denotes the gravitational tensor, and all suffixes not other- 
wise defined denote covariant or contravariant differentiations. 


It is usual to impose the conventional divergence-condition ¢” = 0; 
but this we shall replace by the condition 

¢” = constant. (1.4) 

At a point of space-time where there are no electrons, J, = 0, and 
equations (1.3) reduce, in virtue of (1.4), to 

9°? b.4p+ Gib. = 9. (1.5) 

We propose to solve the partial differential equations (1.4) and 

(1.5), expressing the solution in tensor form, for the case of an 

electron in a space-time of constant curvature K. For this purpose 

we ignore the distortion of the metric due to the gravitational effect 


* The term ‘electron’ is used in this paper in a purely conventional sense, 
and is defined to mean an ideal point charge. 
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of the electron itself, this effect being in general small.* By so doing 
we are enabled to regard the components of the fundamental tensor 
Ju» @8 given functions of the co-ordinates, so that we have a deter- 
minate set of partial differential equations to solve. 

The solution of these equations for an electron in a flat space-time 
(i.e. a space-time in which the Riemann-Christoffel Tensor is every- 
where zero) has already been found,f the result being enunciated as 
follows: 

‘Let an electron e be moving in any manner, and let 7 be its 
proper-time at the point where its world-line intersects the null-cone 
of any point P. Then the electromagnetic potential-vector at P, due 
to the electron e, is l 


é 
bp = Popa (1.6) 


where Cy = Ty/F™ Tims 
and where the subscripts p, q, 7, 1, m on the right-hand side represent 
covariant differentiations.’ 

It is shown that, when space-time is Galilean, the formulae (1.6) 
give the classical formulae for the scalar and vector-potential of a 
moving electron. By an application to special cases, it can, however, 
be shown that they do not give a solution of the equations (1.4) and 
(1.5) when the Riemann-Christoffel Tensor is not everywhere zero. 

The formula which we seek to establish for a space-time of con- 
stant curvature may be regarded as a generalization of (1.6); but, as 
we shall show, it is in fact much simpler, and gives as a special case 
a new formula for the potential of an electron in a flat space-time, 
which is, however, directly transformable into (1.6). 

3efore proceeding to establish the theorem (§ 3) which forms the 
main subject-matter of this paper, it is necessary to prove a theorem 
in the Tensor Calculus. 

§ 2. Theorem 

Let (x®, a1, a, ..., 2”-1) be generalized co-ordinates in an n-dimen- 
sional space of constant curvature K whose metric is given by 
, ds? = g,,,dxtda”, (u,v =0,1,...,m—1). 
be the Riemann-Christoffel Tensor, we shall have 

Bu, ov — K(GyoIpy—Ipv I po)» 


* E. T. Whittaker, Proc. Roy. Soc., A, 116 (1927), 720, § 1. 
+ Ibid. 120 (1928), 1 


Then if B,, 


p> OV 
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and hence also : 3 
a ae 5) 
Guy = 9? Buy, ov a —(n—1)Kg,,. (2.1) 
Let (£°, #1, ..., "-1) be a fixed point of the space, and (x®, 21, ..., a”-1) 
a variable point. Let s denote the length of the arc of the geodesic 
(assumed unique) joining these points, so that s is a function of the 2*. 
Let Q denote the scalar 
Q = cos(K3s). 
Then we shall prove that 
Q satisfies the partial differential equations 
; . ¢ 
Q,, Q« = K(1—Q?), (2.3) 
and Quy cm —_ Q, (nx, vy = 0,1, 2 (2.4) 
where the suffixes on the left-hand side denote covariant or contra- 
variant differentiations, and where (as always) repeated suffixes 
imply summations. 
The proof of (2.3) is simple if we quote the well-known theorem 
that the arc s of a geodesic satisfies the partial differential equation 
> © 95 
$5 = 1. (2.5) 
For since :>= K-arccos Q, 
we have s, = —K-+.(1—Q?)?.Q,, 
with a similar formula for s*. 
Substitution in (2.5) immediately gives (2.3). 
For the proof of (2.4) we take as the canonical form of the metric 
of a space of constant curvature that defined by 


ds? = — K-1(x°)-*[ (dx°)?+ (da1)?-L-. . .+-(da"—1)?]. (2.6) 


By a real or imaginary transformation, the metrics of all such spaces 
can be reduced to this form. 
Suspending temporarily the summation-convention, we have 


Iup = — l K (x°)?, Inv = 0, 


gt! = — K(x°)?, gh®¥ = 0, (u~v;p,v—), 
The only non-vanishing Christoffel symbols are 
{00, 0} = —1/2°; {up, 0} = 1/2*; {Ou, p} = {u0, p} = —1/2°, 
n—1, but p+ 0). 
Substituting in the equations of the geodesics, namely, 


’ da® dab - 
xB, 1} ds ds 
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(the summation-convention being now restored), and integrating, we 
get ultimately the formula 


s = K-are cos] ()2-4 (0/24. "5" (an — any] [2x] 
pal 


for the geodesic distance between the points (x‘) and (#*). 
Hence for this metric, 
1 = = 14\9 
[t is now a matter of straightforward differentiation and substitution 
to prove that 
eg eQ . 
= — a ° — ms 
Oar Gar” ty, a} oa™ Iu 
for 
And since this is a relation between tensors, it must be true for 
all spaces derivable from (2.6) by point-transformations, i.e. for all 
spaces of constant curvature K. 


§ 3. Solution of the partial differential equations of 
electromagnetism 
Consider now a space-time of constant curvature K whose metric 
is specified by 
ds* = g,,datdx’, (u, v = 0, 1, 2, 3), (3.1) 
x° being the ‘time’ co-ordinate. 
Let the null-cone of the point (x) cut the world-line of an electron 
e in the point (#‘); so that, in the language of the older physics, the 
electron was at the point (#!, %, #) at the time Z° when radiation 
left it to arrive at the point (a1, x?, x) at time 2°. 
Suppose the motion of the electron to be given. Then we shall 
know #1, #*, #* as functions of the ‘time’ Z°, say 
a =f(®),  (r=1,2,3). (3.2) 
Let 7 be the proper-time of the electron at the point (%'); then since 
dr* = g,,, ditdz”, 
J, being the value of g,,, at ('), we have, on substituting from (3.2), 
dr? = (dz°)? x a function of <°. 
Integration of this equation gives £° as a function of 7, and hence, 
by (3.2), #1, #, # are also given as functions of 7; thus 
#=#(r),  (r=0,1,2,3). (3.3) 
Let s(x’; €') be the geodesic distance between the points (2‘) and (€'). 
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Then, regarding the ¢' as current co-ordinates, the equation of the 
null-cone of (x*) is 

a(x‘; £*) = 0. 
Since this passes through (Z'), we have 

a(a*; Z*) = 0. 
Hence cos{ K?s(x'; #')} = 1, 
that is, ei. (3.4) 
Q being the function defined in § 2. This equation, combined with 
(3.3), expresses 7 as a function of x°, 21, x”, 2. 

It should be noticed (i) that Q is a function of the 2‘ and of the 
z'; (ii) that in virtue of (3.3) it is therefore a function of the x 
and of r. 

We shall now show how to prove the following theorem : 


The potential of the electron e (whose motion is specified) is given by 


$,= -< = [108 (2) |. (3.5) 


where (i) ¢ is a constant; (it) the partial differentiation with respect to 
at treats 7 as a constant, i.e. ignores the fact that r is given as a function 
of the a‘ by the equation (3.4) ; (iii) + is eliminated after the performance 
of the differentiations by means of (3.4). 

Our justification for describing this as the potential of an electron 
is that in the particular case when K = 0 and the co-ordinates are 
Galilean, (3.5) reduces to the classical formula for the potential of 
a moving electron. 

Before indicating the general proof of the theorem, it will perhaps 
be as well to apply the formula to a particular case. 

Consider an electron at rest* in the space-time specified by 

de? = (1-++-px)-*[de*— = (dz?-+dy*+d22)| : 
(9° = t, 2! = 2g, 2° = y, e* = 2), 
pw being a constant. This is a space-time of constant curvature. 

Since the electron is at rest, 
c= 0=>9=— 2: (3.6) 

] 


hence dr? = (1-++-pt)-*[d72— 
= df?, 


* A certain confusion of language arises on account of the admixture of 
the ideas of the older physics with those of the newer. 


(dz?-+-dg?+ dz)| 


= 
2 
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and we may therefore take 


t=r. (3.7) 
For this space [cf. (2.7)], 
Q = M(1+-mx)(1-+ ma) (1+ pa) (1+ pe) 
, —p*{c2(t—t)?—(y—g)?—(2—2)}] 
; ae eee ee 
by (3.6) and (3.7). 
Hence 


and therefore 


$o Ss ) | = —e/ (3.8) 


$= a, | 8 . (3.9) 
C T 
(3.10) 
Now 7 is obtained as a function of (x, y, z) by the equation Q= 1; 
this gives 
x?+y?+2? —_ c2(t—7)?, 
and thus 7 =t—r/e, where r? = x?+-y?+-2?. 
Substituting in (3.8), we get 
oo = —e/r, dh, = ep/c(1+pa), 4, = 0, ¢, = 0. 
It is a matter of straightforward algebra to show that this is indeed 
a solution of equations (1.4) and (1.5) for this particular metric. 
It will be seen that the constant c has the dimensions of a velocity, 
and can (in the case considered) be interpreted as that of light. 


§4. Proof of the theorem 
A complete proof of the theorem of the last paragraph cannot be 
given on account of its length; but the following should serve to 
indicate the method adopted. 
Ignoring the multiplicative constant —e/c, we have to show that 


the vector defined by 
aT, /2@ 
2 log (&& 4.1 
Fu Oa [ °8 ( Or )] ai 


satisfies the partial differential equations 
9? db rapt F6 be = 9, (4.2) 
pi = constant. (4.3) 
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This may be done by direct substitution of (4.1) in the last two 
equations, making use of the results of § 2, namely, 
2,9" = K(1—Q”), (4.4) 
and Quy = —Kg,,9. (4.5) 
A certain confusion may arise on account of the fact that Q is 
a function of both the x‘ and of 7, where 7 itself is a function of the 


x defined by the equation Q=1. (4.6) 


Covariant differentiations of Q with respect to the 2‘ which treat 
7 as a constant will be denoted by the simple addition of suffixes. 
Thus oQ 

a= aes (7 constant). 

Covariant differentiations which treat 7 as a function of the a‘ defined 
by (4.6) will be indicated by the bracketing of the corresponding 
suffix ; thus 
eq 


a (7 not constant). 
ox 


Qu) as 


Hence Q,=@ 4.20 a 
a 

— Q.+9,7, (4.7) 

where the suffix 7 denotes partial differentiation with respect to r. 
With this notation we have to show that the vector ¢, defined by 
Pu —- Qu Q,, (4.8) 
where Q_, = a "@ —— at 
" + OtPOr dOrox 
9°PD (0x) (B) + Gs $. = 


and gd u(4) = constant. 


. , Satisfies the equations 


Equations (4.4), (4.5), (4.6) remain unaltered. 
Differentiating (4.4) twice with respect to 7, we get 
Q.9" — KQQ,, (4.11) 
Qe. Ql + Qn! = —KQQ,,—KQ?. (4.12) 
In (4.11) use is made of the fact that Q,Q¢=Q+Q,,. Differen- 
tiating (4.5) with respect to 7, we have 


Quy — —Kg,,9,. (4.13) 











THE POTENTIAL OF AN ELECTRON 
Also, differentiating (4.6) with respect to x*, we have 
Qn) = 0, 
and hence, by (4.7), T, = —Q,/Q,- 
Finally, on raising the suffix v in (4.13) and contracting, 
Qt, = —4KQ,. 
Now Ody. 
iNOW Pui ~ mae Onat or Tx 


— (Q54/Qr)at+(Qry/Qr)2T a» by (4.8), 
= QF (Qruat C.utd— Q77(QQrat Qiu Vr77 0): 
On multiplying by g#“, summing for w and «, and using (4.14), 
g” “P u(a) = Q; : — Q; *(Qrr e+ Qu 0%)+ Q7°0,.9"Q,,- 
Using (4.15), (4.12), and (4.11) we at once get 


9" “bu(a) = —3K = constant, 


so that ¢,, thus defined, satisfies (4.10). 

By a similar process, using (2.1) (with n == 4), and remembering 
that we may put Q = 1 after all differentiations have been performed, 
it can be shown that the expression given for ¢,, also satisfies (4.9). 
The work is lengthy, but straightforward if the distinction between 
bracketed and unbracketed suffixes is borne in mind. 


§5. The case K= 0 
We now deduce a formula which is valid for a flat space-time, 
i.e. a space-time of zero curvature. 
Let P be the scalar defined by the equation 
P= (1—Q)/K (5.1) 
= 2{sin?®(}K*s)}/K. 
Then lim P = }s? 
K—0 
, =Q, (5.2) 
where Q denotes one-half the square of the geodesic distance (in a flat 
space) between the points (x*) and (z*). 
By (5.1), Q=1-—KP, (5.3) 
and substitution in (4.4) and (4.5) shows that P satisfies the equations 
P, Pt = 2P—KP*, (5.4) 


P., =9yA1—KP). (5.5) 
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Making K +0, and assuming that 


ae of. 
lim — = = lim P }, 
K-00x¥ = 0x" \ K-+0 
we find that the scalar Q as above defined satisfies the partial 
differential equations 0,0 = 20, (5.6) 
Quay = GIpr (5.7) 
Furthermore, (4.6) becomes 
Q=— 0. (5.8) 


€ 


Also, substituting from (5.3) in (3.5), and proceeding to the limit, 


$, = lim[—eP,,,/cP,] 


K-—0 


A 


in the notation of § 4. 
Hence, for a flat space, 


Fu 


which may be written 


$, — ii [ton = (5.9) 
With the help of (5.6), (5.7), (5.8), this formula can be transformed 
directly into (1.6). It has the advantage of being considerably 
simpler than (1.6), since its application involves only simple partial 
differentiations instead of an elaborate series of covariant differentia- 
tions. 

Although we have deduced (5.9) by a limiting process from (3.5), 
it can nevertheless be established independently by a method similar 
to that used in proving (3.5). 

Application. 

Apply (5.9) to the case in which the co-ordinates are Galilean, that 
is, when the metric is given by 


ds? = dt®— ~ (da? 4-dy?-+-d2?). 
5 


Then, of course, 


* These equations are established independently in a paper to appear in 
Proc. London Math. Soc. 
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Hence 
aQ = eQdt aQ dx eQ dij @Q dz 
Or Ohdr| d& dr | dy dr & dr 
1 - ™ i 
= —(t—fey+ fe), + Yat (@—2)r3} 
= ys, say, 
dt 
PA > 
Hence, by (5.9), 


where Vo = 


eV 


oh (5.10) 


ed 
=-—-— lo ) = 
bo - na gy 
with similar expressions for ¢,, $2, $3. 
Now 7 is given as a function of 2, y, z by the equation Q= 0, 
which gives the familiar relation 
where 


Hence ‘= 


by = —e/[7—{e—B)0, + Y—Dry +204], 


with similar expressions for ¢,, $9, ¢3. This is the classical formula 
for the potential of an electron which was moving with the velocity 
(v,, Vy, V;) at the point (%, 7, Z) at time ?, the radiation then emitted 
arriving at (x, y, z) at time f. 





QUADRATIC EQUATIONS WITH AN INFINITY 
OF UNKNOWNS 
By U. 8. HASLAM-JONES AND E. C. TITCHMARSH 

[Received 20 May 1930] 

1. THE theory of linear equations with an infinity of unknowns is 

a standard branch of analysis, and the subject of a well-known book 

by F. Riesz. Equations of degree higher than the first seem hardly 

to have been considered in detail. They are no doubt less important 

and more difficult to manage than linear ones; but there are some 

cases, at any rate, which, as we shall show, can be dealt with by 

quite simple methods. 

In this paper we attempt to find conditions for the existence of 
real solutions of an infinity of quadratic equations in an infinity of 
unknowns. The most general form of the equations which we con- 
sider may be written as 


> > Annvémén=Cy Ww=1,2,3,...), (1) 
1 


m=1n 
where the nuinbers A,, ,,, and C, are all real. 

The problem, so stated, is formally analogous to the corresponding 
problem involving linear equations; but many of the results of the 
theory of the latter depend upon methods which are available for 
the linear case only, and do not extend to the solution of quadratics. 
For our present purpose the method of successive approximations 
appears to be the most profitable. 

2. We assume that the numbers A,, , and C, are all positive. 
Then by means of the substitution 


b= (4.".)* 


and division by C,, the equations (1) are reduced to the standard form 


« 
> Qn, n,v& men = 1 (v= 1,: 
m=1n=1 
where a, , , = 0 for all values of v. 
> 


We take, as a first approximation to the solution, x7 ,=—1. The 
successive approximations are then defined by the formula 


F + > z Amn, n, v&m,r—-1 © n, r- _— 1 ( =1 2, ? a} (3) 
m nt 


the sign of x, , being the same as that of 2, ,. 
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Our main theorem is as follows: 
THEOREM. Let 0 be the upper bound of the sums 


Pes Men wel (v= 1, 2,3,...), 
mn 


and let any arbitrarily chosen sequence of signs be given to the variables 
x,. Then there is a number a, between 33 and }, such that, if 0< a, 
the equations (2) have one, and only one, vl bounded solution with the 
given signs. The solution is 

x, = limz, , 


ro 
In the first place we restrict ourselves to solutions in which all the 
variables are positive, and so assume that every 2, , is of positive sign. 
The proof consists of a number of stages. 
If 0<4, the numbers x,, are recl and bounded for all values of 


v and r. 


Let M, = Max(2z? ,), 


so that, in the first place, M@, = 1. Then from (3) 
Ly, tyr S14 2D lem, n, ol Sn, 1%n, r- Bee <14H,_ 12> [Om, nv 


mn 
and so M,<1+6M,_,. 
Starting with M, = 1, repeated application of this inequality shows 
that l 
s, <14+6+...4+0-1< <j 9 (4) 
for all values of r. 


Hence also 
Zz B+ F 3 > ln, n,v Xm, r-1 vn r-1 


> 1—0M,_, 
~~ 6  1—20 
1—@ 1-86 
by (4). Combining (4) and (5), we have, if @< }, 
1—20 1 
x ‘ 6 
i < v,r <, - | ( b) 
It follows that the numbers 2, , are real and bounded. 
4. Let « be the number such that, if 0<«, x, , is real and tends to 
a limit X, for every value of v, while if 6 > « this is not necessarily the 
case. Then «x lies between }3 and 4 
3695 M 


0< 
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From (3) 


”» 


[ee +a, | > lA a, vite, 2 Un, 9 See, v2 Xn, r-1)| 
n 


< 0Max(|z,, ,-%, -— 


"m7 "n, f m,r—-1 Ln, r—1|)- 
m,n 


We write Max(|2; -—2?, 
; 
and we have, if m #«n, 


Xm,r Ln, r —ZLm, r-1 Xn, r—1| 


— 5 (Tm, a Xm, r-1(Lp, rp, r- 1) Si (Xp, r+2p, r-1)(m, r— Xm, yg) 


| Xm, rT a 


a . | . 

m,r—1 (2 3 1 nr Tn, r-1 1.9 2 

2 j (x, y “~~ r- i)+ ( ee 
1Tn,rTUn, ra Os a 


. | ” . | . 
5 (Gnetinet 4 Saree ‘) 
Ln, rT&p, r-1 Xm, eta m,r—1 
expression in brackets is of the form 
xr+-, 
x 
and so has its maximum at one end of the interval over which x 
ranges. Since by (6) 
] 
/ ¢ - - 
4(1—20)<24<— : 
( V(1— 26) 
it follows that 


_ | 
Um, rn, r—&m, r—1 © n, r-1| S 7 et ~ 


If m =n we have of course simply 


9 9 1 
|x ’ —Tn,r-1| < ) 


& 
Therefore from (7) 


Wi—20)4+—_" | —)8., say. 


Br < 101 ; 
~ | V(1—26)} 


If A< 1, the series & 5, converges, and so every series 
2 2 
(%), -—3 »,r—1) 
r=2 
converges; and hence 2,, tends to a limit, X, say, as r tends to 
infinity. Also from (6) 


1— 20 xc! 
19 “"”" “1—9 
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The sufficiency of the condition A < 1 shows that the number a is 
at least as great as the positive root of the equation 
, ‘Re 
yo! (1—20) + , i. 8 
40 | (1-20) + a3 (8) 
This root is slightly greater than }5, since, if we put 0=}5, the 
left-hand side of (8) is equal to 33, and so is just less than 1. 
On the other hand, a simple example shows that a cannot be 
greater than }. Taking the system 
a?+Oe3=1, 223—O22=1,... 22?—Ox2_,=1,... (9) 
we have, for v>1, 


a? = 140408+...46'-1, 


ry, 
while af, = 1—62z3 4 = ‘1—6—@—...—6r-1, 
Thus, if 6>4, a7, is negative if r is sufficiently large, and the 
approximations cease to be real. 
” 5. It is easily shown that the numbers X, are actually solutions 
of the equations (2). In fact, more generally, if 6 <1 and z, , is real 
and tends to X,, as r tends to infinity, for every v, then the set of numbers 
(X,,) ts a solution of the equations (2). 

By our initial assumption, the series 

T¥ la | 


a im, nv 
mn 


is convergent for every v; and, by (4), the numbers x, , are bounded, 
since <1. Hence each of the series 
3's 


m 


a a ae 


m,n,v“m,r 


a 
7 


converges uniformly with respect to r. Therefore 
> ee "> = : ; 
a An, n, - = Xi, tna lim an Amn, n, yp&m, r% n,r 
; 2 . 
= lim (l—25, ,,3) = 1—X}, 


mn ren im n 
rr 


which is the required result. 

6. We have thus shown that, if @ is less than the positive root of 
equation (8), the method of successive approximations leads to a 
solution (X,) of which all the members X, are positive. In order to 
find a solution in which the members have any arbitrarily chosen 
sequence of signs, say x, has the sign of (—1)*”, we begin by making 
the substitution x, = (—1)*2’. 

M2 
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The problem is thus reduced to that of finding a positive solution 


of the equations LTT a epee) 
a v + > ea Amn, n, yt mt a 1, 
mn 


where |a a Since 6 is unchanged thereby, a solution 


’ ' 
m,n,v m,n,vi* 


with any given signs exists under the same conditions as the solu- 
tion of positive sign. 

7. The general question of the uniqueness of the set of solutions 
so obtained seems to be difficult to answer. In the case of linear 
equations solutions can, under certain conditions, be found containing 
an arbitrary parameter. Naturally, corresponding systems of quad- 
ratic equations exist. For example, if 


xi—OaZ=1,... 25—62},,=1, ... (10) 


we can give 2x, any value, and then obtain the values of 29, x3, ... in 
succession. Hence the most general solution can be put in the form 


cok. . 2 
- 110° 


a 
where € is arbitrary. Thus, however small @ is, the solution of our 
system of equations is not necessarily unique. But, if 0<@ <1, the 
solution is bounded as v-oo only if € = 0; and then the solution is 
that which would be obtained by the method of successive approxi- 
mations. 

We can show generally that if 6 is less than the positive root of 
equation (8), the solutions obtained by the method of successive approxi- 
mations are the only bounded solutions. 

Let (Y,) be a bounded solution of the equations (2), say |Y, 


Y,? < i-}- > p 2 Gn, n, v | ) < 1+6M?. 


— 


<M. 

Then 
mm 

Starting from this new inequality, and repeating the process, we 

obtain 


2 < 1+0(1+0M?), 


and generally 


Using this inequality, we have 
7 | > > wz ,~ 1-20 
LY, *>1—2> a ee St 
mn l —— 
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The proof now follows the steps of § 4 above, for, as shown in § 6; 
we may assume that every Y, is positive. We have 


<|>>da Y.. ¥.1< 


m,n,v~m ? 
mn ™ San 


—*X 26 
| < 
Cy, r-1 eer i 
we have, as in § 4, 
| 72__ 42 | sy > a 4 | 
|} v I< 22, (9 1m, n, AYnYa— Xm, r— -1%n, r-1)| 
mn 


Vv, 


r-2 
r-16 
1-8" 
Hence, by induction, (11) is true for all values of r. Hence, if A <1, 
, tends to Y, for all values of v as r tends to infinity. That is, 
X,=Y, 
and Y, is identical with one of the solutions already obtained. The 
proof of our main theorem is now complete. 

8. It remains to decide how nearly the conditions obtained are the 
best possible of their kind. Simple examples show that, if 6 > 4, not 
only may the method of successive approximations fail, but also real, 
bounded solutions of the equations may not exist. Thus the most 
general solution of the system (9) above is 

4 
oo (—>+ ef <a e } 


( = 
‘ke 2)+ am} 


ae |e me || 
where é is arbitrary. This solution is real and bounded if, and only 
if, €=0and 0< }. 

The question may be stated in a rather different way. We have 
already defined the number a to be such that, if @<«, the method 
of successive approximations succeeds, while if 6 >a it may fail. 
Let 8 be the number such that, if @<, the equations have a real 
solution with any given system of signs, while, if 6 >, this is not 
necessarily the case. It is then clear from the definitions that B > « 
but we do not know any lower limit for 8 apart from that for «. On 
the other hand, it follows from the above example that 8 < }. It is 


of course quite possible that « = fB = }. 
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9. There are similar problems in the case of » quadratics in n 
variables, where n = 2, 3,.... In each case we can define numbers 
n> Bn: in the same way as the above a and f. Thus a, is the number 
such that the approximation method succeeds for the equations 


x 


x?+ary+by? = 1, 
ca?+-daey+y? = 1, 
if 6 = Max(|a|+ |b], |c|+|d|) <a, 
but not necessarily if @ > «,; and similarly in other cases. 
Since the system with m unknowns can be considered as a particular 
case of that with n--1, it is clear that 


he Oe nS Oye nse 
and that Bo Ph SE, Pee ee 


and of course f,, > «, for every n. 

The case n = 2 is of some interest because it is quite easy to deter- 
mine f, exactly, and because here a, is actually different from fy. 
We find that 8, = 1, while a, << }(V5—1). This last fact is seen by 





considering the equations 
xi+623 = 1, x3— Oa = 
The approximations are x, , = 1, x, , = 1, and 
ay 9 = V(1—6), Ly 9 = V(1+8), 

% 3 = V(1-—-0—6), Xp 3 = V(1+0—6?) ; 
and 2, , is imaginary if 1—@—6? < 0, ie. if @ > 4(v5—1). 
To show that 8, < 1, we observe that the solutions of 

x?—Oa3 = 1, 22—Oxi = 1 
are complex if 6 > 1. 

10. The statement that B, >1 is equivalent to the following 

theorem. 

If \a\+\b| <1 and |\c|+|d| <1, all the solutions of the equations 
x*+ary+by? = 1,) (12) 
cx?+-day+y? = 1 ) 

are real. 

We give two proofs. 

(i) Writing y = xr, we have 

x2(1+-ar+br?) = 22(c+dr+r?) = 1. (13) 
Hence, since x = 0 is not a possible solution, 
ltar+br? = c+dr+r?, (14) 
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and, since 1—ec and 1—6 are both positive, this equation has two 
real roots. 

If p is a root of (14), 1+ap+bp? and c+dp-+p? are both positive. 
For suppose on the contrary that 

1-+-ap+bp? < 0. (15) 
Then |p| > 1; for otherwise |ap+bp?| < |a|+|b| <1, and (15) is not 
satisfied. Therefore 
c+dp+p? > p*(1—|e|—|d|) > 0, 
which is inconsistent with (14) and (15). 

It follows that (13) is satisfied by two real positive values of 2°, 
and so the equations (12) give four real pairs of values of x and y. 

(ii) Considering (12) as the equations of two conics, any conic 
through their intersections is 

x*+-axy+by?—1-+-k(cx?+day+y*—1) = 0. 

This is a pair of straight lines if k = —1, or if 

4(1+-ke)(b+k) = (a+kd)?. 
This gives real values of k if 

(1—bc)? > (d—ac)(a—bd), 
and the straight lines are real if either of the inequalities 

1+ke 0 b+k 

i+: * 1+k 
is satisfied. These therefore are the conditions that the four solutions 
should be real. Now 

0 <1—|a|—|b| <1—|a|—|6|(\e|+-|d}), 

so that |ja|+ |bd| < 1—|be}, ja—bd| < 1—be. 
Similarly, |d—ac| < 1—bc, so that (16) is satisfied. Also if |k| <1 
the first inequality in (17) is true, while if |k| > 1 the second is true. 
Hence the result. 

11. The case n = 3 is more like the case of an infinity of equations. 
Consider, for example, the system 


2 9, 


u 


aj +6a3 = 1, x3—0a3 = 1, 22—Oar3 = 1. 
These give oo. 1 P 1 
a i7* * oe 
and for real solutions we must have 6<}. Thus f; <4. We have 
not, however, determined the actual value of f;, or of any of the 


other constants. 








ON BOUNDED BILINEAR FORMS IN AN INFINITE 
NUMBER OF VARIABLES 


By J. E. LITTLEWOOD 


[Received 10 April 1930] 


1. Proressor P. J. DANTELL recently asked me if I could find an 
example of a function of two variables, of bounded variation accord- 
ing to a certain definition of Fréchet, but not according to the usual 
definition. 

The problem is equivalent to the following: to find a bilinear form 


, a a . 
Q(x, y) “a Zz a Ann* mYn (1) 
in an infinite number of variables x, y, such that > > |a,,,,| is diver- 
gent, but Q is bounded for all x, y belonging to the range S defined by 


x.|< 1 (= 1,2, ...), YnlS 
We shall say that ‘Q is bounded by H in S’ if 
M N 
. 2 Oot eee 
m=1 n= 
for all x, y of S and all M, N.+ 
There is a well-known theoryt of forms Q(u,v) bounded in a 
‘Hilbertian’ space 
> Mnl?<K, > |)? <K, (2) 


and our problem is answered, in the affirmative, by Hilbert’s form§ 
1 UnUp 
3 2 m—n’ 
which is bounded for wu, v satisfying (2). Calling this Qo(u, v), let us 
fix a set of u, v satisfying (2), and consider 


U(x, y) — Qo( UL, vy) = > > Ginn Xm Yn ? 


where 
U,V 
Ann = m * (m 
m—n 


Evidently Q is bounded in S. On the other hand > > |{a,,,,| is 
divergent if, for example, we take 


=n), == © (on = #). 


7 


Um = Vp, = m-*(log m)-%, L<a<l. 


+ We shall see in § 3 below that the double series is then convergent. 

{t For an account®of this, see F. Riesz, Les systémes d’équations linéaires 
a une infinité @incornues (1913, Borel collection). I shall refer to this book 
as ‘Riesz’. 

§ The dash denotes that the terms for which m = n are omitted in the 
summation. 
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The form thus resulting has > > |a,,,,|\!** convergent, and it is 
natural to inquire how far removed from convergence the series 
> > |@n»,| can be for a bounded form. A complete answer to this 
problem and a proof that it is complete involve solving also the 
converse problem, that of necessary conditions for a form Q to be 
bounded in 8S. Both problems are solved, in a sense completely, in 
the present paper, in writing which I have profited by a discussion 
with Professor Daniell. The converse problem, of necessary condi- 
tions, turns out to be the more interesting of the two. It is inci- 
dentally quite unlike the corresponding problem for Hilbertian space, 
where the results lie more on the surface.t 


Let 2 AY @ . 
b,, aac >. Ann | ? Cm _ a IGnn "ys 
m=1 n=1 


i 4) 
2 em = C, (> > | 


it)? = D. 


a mni 


> b,, = B, 
n=1 m 
(B, C, D may of course be +00). Then I prove 

THEOREM 1. (1) Jn order that |\Q|< H in S it is necessary that 


B, C, D should be less than AH, where A is an absolute constant. 
(2) Given any positive Ps Ins Tmn for which 


lim p,,= limg, = lim r 


m2 nn m,n->2 


there exist bounded Q for which 
- : < ae i 
p 3 Pp b > YmEm > — Tmn [Amn |* 


n° n? 


= 0, 


mn 


are divergent. 
2. Preliminary remarks. (i) Since we are free to take 
UM4i +: =r%y => 0 
(if e.g. N > M) we may replace M by N throughout the definition 
of boundedness. 
(ii) For fixed M,N |Q),,n | attains its upper bound in S. Since 
M 
Qu = 2 Xm Yn» 


m= 
where 


N 
Y m _ Zz Ann Yn: 
n=1 


+ They are that the series > |a,,,| and & |a,,,|* converge to bounded sums. 
m 


n 
They may be deduced at once from the results for space S, but the converse 
deduction is not possible. 








166 J. E. LITTLEWOOD 


M 
Qin and |Q),,x| attain, for fixed y’s, the (same) maximum }> |Y,,| 
1 


when z,, = sgn Y,,. When |Qj,x| attains its upper bound all relevant 


M 
lz! and |y| have the value 1, and > |¥,,| and the corresponding 
N : 
expression > |X,,| are each equal to the upper bound. 
1 

(iii) We may distinguish three forms of our theorem in which 
respectively (1) the @,,,,; % > Y, are complex, (2) the a,,,, X,, Y, are 
real, (3) the a,,,, are real and z,,, y,, are each restricted to the values 
+1 and —1. After remark (ii) it is easy to see that the three forms 
are essentially equivalent, and we shall confine ourselves for the most 
part, and unless the contrary is stated, to case (1). 

(iv) Our results are most easily understood when stated, as in 
Theorem 1, in terms of the ‘convergence’ of the series B, C, D; but 
their true character is best brought out when we consider truncated 
forms and series in which 2,,, y,, are zero (or, alternatively, a,,,,, is 
zero) when m > N orn>WN. We prove, e.g., on the one hand that 
if |Ovyv|<H in S, then > > |a,,,,|'< AH, and on the other that 


m,naN 
there exist forms Qyy (with arbitrarily large N) such that 
\Qxwl <A (> > lun!) 
m,naN ] 
in S.¢ It is in fact this type of result that we shall be concerned 
with up to the very last moment. 

(v) It follows from the familiar identity 

Q(x, y) = OF(@+-y), 3(x+y)}—Q[h(a—y), H(x—y)}. 
that the bilinear form Q(a, y) is bounded in S if the quadratic form 
Q(x,x) is; and the converse is trivial. Our results have therefore 
obvious applications to the theory of bounded quadratic forms, which 
need not be gone into in detail. 

3. Whether the double series (1) converges in S has really nothing 
to do with the question of boundedness. But since the results about 
convergence are clear cut and help the meaning of boundedness, it is 
perhaps desirable to give them, though the proofs are a little tedious. 

If Q is bounded by H in S the double series (1) converges uniformly 


+ It is possible further to have all the |a,,,,| equal. 


t For the corresponding problems in Hilbertian space (where the arguments 
are more interesting), see Riesz, 81 et seq. 
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in S in Pringsheim’s sense (naturally to a sum Q satisfying |\Q| < H). 
Conversely, if the series (1) converges everywhere in S, the form Q is 
bounded (so that the convergence is necessarily uniform). 
The first part may be seen as follows (the argument is best followed 
with the help of a diagram). The series > |a,,,,| and > |a,,,,| must 
n 


m 
evidently converge for each n, m respectively. Hence, given any rv, 
there exists a A(v) such that 
> > lAnn| <€, > > l4nn| <€. (3) 


n<y m>X mv n>X 
Let H be the upper bound of 
mn 
Qin |= I> > Amn Xm Yn | 
for varying x, y in S and varying m, n. Then there exists an L such 
that Q,;, exceeds H—e for suitable x,y in S. Let A=A(L), M>A, 
N >A, and let us take any set (x, y) in S for which 
™=0, y=0 (L<I<A). 
M -N 
Qun=@rr+ > 2X +R, 
m=A N=A 
where, in virtue of (3) with v= A, 
|R| < > > l@nn|+ > 2 lAmn| < 2e. 


n<Lm>a m<L n> 


Then 


The variables x, y in the first term in the right-hand side of (4), and 
those in the second term, are completely independent. Let us 
choose the first set so that Q,;, > H—e,t the second set arbitrarily. 
The second term cannot exceed 

H—(H—e)—R < 3e 
algebraically. Since this is true for arbitrary variables (in S), we 


must have 

M N 
> ¢ AnnXm Yn < 3e (Mu, N > A) (5) 
m=A n=A 


for all x, yin S. If now A, = A(A) we have in S 
Ni 


— = : we ' + 
~ Ann®m Yn S > > Ann | <€ (N 1> A,), 
Ay 


— — 
A n=: MSA NPAT 
M, 
» 3 Ss Fg © 
2 2 Wn *m9Yn S 2 a 
nsaAim=Ay nsA mPpaAy 


Combining (5) and (6) we have 
\Qar.ni— Qari < 5e (M,, N, > A,) 
+ Ifa bilinear form takes a value qg in S it takes also the value —q; we 
have only to reverse the signs of, say, all the y’s. 


l4nn| <e (M, Pe A,). 
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in S, and this is a sufficient condition for uniform convergence in 
Pringsheim’s sense. 
It follows incidentally that given ¢« there exists a v = v(e) such that 
Q—Q,, is bounded by « in S; the form Q is ‘completely continuous’ .+ 
The converse is easier. If (1) converges for every set x, y of S, 
the rows and columns must evidently converge absolutely; let G, 
denote the sum > > |a,,,,| taken over the gnomon of pairs (m, n) for 


which either m <v orn<v. If now Q is unbounded there exists an 
N, and a set (x, y) in S with suffixes not exceeding N,, such that 
Qn, Nn, (e, y) > 1; then an N, and a set (x, y) in S with suffixes ex- 
ceeding N, and not exceeding N,, such that Qy, v,(2, y”) > 14+ Gy,; 
then an N, and a set (x®,y®) with suffixes exceeding N, and not 
exceeding N,, such that Qy, v,(e®, y) > 1+G@y,; and soon. Let now 
(a*, y*) be the set in S obtained by combining (2, y™), (a, y®), .... 
Evidently 
Qn, n,(z*, y*)—OQn,_, ny-y(4*, 9*) > Qn, w le, yy”) —Gy,_, > 1 

and (1) does not converge at («*, y*), a contradiction. 

4. I come now to our main problems, and begin with the first part 
of Theorem 1. The inequality D < AH is the most interesting of the 
three necessary conditions for boundedness; but it is a consequence 
of the other two and the ‘pure’ inequalityt 

D< A(B+C), (7) 
and it is convenient to begin by proving (7). In this we may suppose 
that the a,,,, are real and non-negative, and that B, C are finite. If 
we interchange a pair of rows n = 7, n, of the matrix (a,,,,) we do 
not alter the value of any c,,, the new b,,,, 6, are the old ones inter- 
changed, and the other b,, are unaltered. We can therefore permute 
the rows-so that the c,, are unaltered and the new 6b, are the old 
ones rearranged in descending order. We can now permute the 
columns so that the b, are unaltered and the c,, are rearranged in 
descending order. Since these permutations do not alter B, C, D, it 
is enough to prove (7) in the new case in which c¢,,, 6, are non- 
increasing. In this case we have 

nb, < B, Mem <C. (8) 

+ Cp. Riesz, 94 et seq. 

t Iuse A throughout to denote a positive absolute constant, not necessarily 
the same from one occurrence to another. 

§ This is possible even when there are an infinity of 6,, because B is finite 


and b, — 0. 
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3y Hélder’s inequality 
Sahn <( 1) E Cahn!) <n10h < BY, 
man m< 1 


man Gi 


x 2 Gan < Bt. 


nmon 
Similarly 
— > Sat,=> > <et. 


Hence Dt < Bt+ct 
D<A(B+0C). 


5. Our essential problem, then, is to prove the inequalities B < AH, 
C < AH. It is enough to prove these for truncated forms Qyy and 
an arbitrary N (when 6?, B,c2, C all become sums of a finite number 
N of terms).¢ The proof depends on a ‘mean value’ principle (which 
I have used before)t, and on a special application of it that seems of 
independent interest. I set out the argument in Lemmas | to 4, 
of which Lemmas 1 and 2 are well known. 

Let A> 0 and let M)(f) denote the Ath root of the mean value of 
the number |f| with respect to a set of discrete or continuous 
variables. The means we are principally concerned with are taken 
with respect to the corners of the ‘cube’ in NV dimensions, 


—l< Xm <1 (m < N), (9) 
in which case uX(f) = 2-N F Iflas, %,..- 520) P 
the sum on the right being taken over all distinct sets (21, % ,...,%y) 
in which every x is +1 or —1. 


Lemma 1. M)(f)<M,(f) (0<A< 
Lemma 2. LetO<a<B<y. Then 
M,(f) < M3(f) M}-*(f), 


where 3 __ay—B) ¢ 
wo —a) . 


Bly 


+ For Qyy is bounded by H in S if Q is, and e.g. B(Q) = lim B(Qyy). 
; yn 
t J. E. Littlewood, ‘On the mean values of power series’, Proc. London 
Math. Soc. (2) 25 (1926), 328-37. 
§ Lemma 1 is a case of Hélder’s inequality with index y/A. Lemma 2 
».g. for means over a discrete set of N terms) is a consequence of 


1 , 
ME(f) Xl F|"|F|F- u < (SEs) cs] evelyn 
with up = a, oo age =¥. 
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LemMA 3. Suppose that M,(f)>T, M,(f}<AT.+ Then 

M,(f)< AT, M,(f)>AT. 

The first inequality follows from M, < M, (Lemma 1). The second, 
an immediate consequence of Lemma 2, is our “mean-value principle’. 

Lemma 4. Let M’s denote mean values with respect to either (i) the 
corners of the N-dimensional cube (9), or (ii) all the points of the 
N-dimensional cube (9). Thent 


M (04% + a9%_+-...+- ayy) 


lies between two constants A. 

In geometrical language: the mean perpendicular from the corners, or 
from all the points, of the cube (9) upon any [N—1] (> ¢,,%, = 9) 
through the centre lies between two constants A, whatever the value of 
N and whatever the position of the [N—1}.§ 

We may suppose > |a,,|?= 1. Then 


Mi > af.) = >, 0, ,0,0,\ = > a,0,M ,(2,2,) 


mm rs 


(10) 


since M,(x,x,)=0 if r~s and M,(x2) is 1 or } according as the 
means are a type (7) or (2). 

Next, 
Mi > Xn) * M,/ 


rn eo ee eS 
Mp4, My) = 2. XO 5,0, IM, (21,0.0,20,,). 
r,8,t,u 


r, 8, t, u 
Now, for fixed r, s (equal or unequal), M,(x,%,x,x,) is zero unless 
t=r, u=s, or else t= 8, u=—r. Paiepibigte 3 
“ys » )<2> Jaa IT Int? ¢ 
Mi(> XmXm) <2> - -Xe\|= . (11) 
r,smiNn ", 8 
Lemma 4 now follows from (10), (11), and Lemma 3. 
We may observe in passing that a similar argument proves that 
for A Be 0 : *~ Z 
My( Z, mt a 
+ 2 
‘ ( — Xm ) 
lies between two positive constants A, depending only on A. In fact, 
+ In our application the A here is, as it should be, an absolute constant. 
The «,, may be complex. 
§ It seems an interesting problem to determine, for given N, the positions 
of the [V—1] that make M, a maximum and a minimum. 











; ON BOUNDED BILINEAR FORMS 171 
if 2k is the even ae next above A-+-1 we can prove, much as 
before and supposing > |a,,|? = 1, that 

Myla OpiXm) < Ay = A). (12) 
Hence M) < My,,< A) by Lemma 1. For the opposite inequality we 
have M, >M,=A if A >2; while (12) and Lemma 3 with «=A, 
B= 2, y=4 give M, >A) if 0<A<2. 

6. Consider now part (1) of Theorem 1: it is enough to prove 
B < AH for a Qyyn bounded by H in S. Let M’s denote means with 
respect to the corners of the cube (9). If (2, 2, ..., y) is a corner 
(so that each 2 is +1 or —1), then for a suitable set of (complex) 


yin S a N 
\Q| = = XW=f, 


m=1 


say, where X, => a Hence f < H for all corners. 


m=1 mum" 
Now the mean of a sum is the sum of the means; therefore 
H > Mf)=> My(\X,)). 
Hence, by Lemma 4, : 


a>z A (2 | Gn") = AB, 


the desired result. 
Turning now to part (2) of Theorem 1, I show first that there 

exist arbitrarily large N and forms Qy yn for which 
’ \Qnn| <AB, \Qnn| < AC, \Qnni <AD (13) 
in S. 

Suppose that (Cnn) (m0 = 1,3,...,2) 
is a set of real direction-cosines of an orthogonal set of axes in 
V dimensions, or, more generally, that the e,,,, are complex and 


satisfy P 
utisly Se ae “lores (14) 


a Crn€ gn O(r a s) s 
Fé 


n 
Let now Q= Qyn; ann = N- €mn (m,n < N). 


Then in S 
aol< SIS N-I m/s -1N 2\4 
(| “SS Ze > N Cmn& m < N { = N a > mn Tn ) ? 
nim n m 
<V ; fot ee Ss 5 
N |Q|? “SS Los Le Cry Xe snXs = a Cpl g a Crn€sn 
nr,s r,s n 
= > 223..=> |e. fad 
r 


rvs rs ri ? 


IQr?< 


r,s 
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The numbers B, C are each 1, so that the first two inequalities in 
(13) are satisfied whatever the system of direction-cosines (e,,,,). It 
is, however, the third inequality, concerning D, that is important. 

> lemn|*)? = N*M (enn) 

(the mean being over m,n < N). Now M,(e,,,,) has the fixed value 
N-?; hence, by Lemma 1, D cannot exceed 1. D actually attains the 
| are all equal; and we naturally try to make them 


Tre Y 
Here D=N-\ 


value 1 if the /e,,,, 
as equal as possible. 

Equality of the |e,,,,| is obtainable, for any value of N, if the e’s 
are allowed to be complex. We may take for example the system 

seid 
enn = He * ? 
or a system obtained by permuting the rows and columns of this one. 
If N-+-1 is a prime p, the system 
Emn = N-*x,(m), 
where x,(m), x.(m), are Dirichlet’s characters to modulus p, is ob- 
tained by such a permutation. 

If the e,,,, are to be real direction-cosines it is not possible to secure 
equality of the |e,,,,| for every value of N (e.g. it is impossible for 
N =: 3). It is, however, a fact (presumably well known, though I have 
not met it before) that in space of N = 2” dimensions it is possible 
to have a set of N orthogonal real vectors (e,,,,) for which |e,,,,| = N-}.t 
In 4 dimensions there is the elegant 

—1 1 1 

1 —1l 1 1 

1 1 —l ] 

1 1 1 —1 
(in which we have suppressed the constant factor N~-*), and we can 
form the matrix M7, , for VN = 2”*1 from M, by juxtaposing matrices in 
the manner vu Py ) 

ve LE =~ 

+ That Q can satisfy the first two is obvious at once from the trivial form 
in which a,, = Nim =n), = O(m 3 n). 
The importance, such as it is, of these inequalities is that they can be satisfied 
at the same time as the third. 

t Mr. H. D. Ursell, to whom I mentioned the problem of making the lennl 
as equal as possible, discovered this independently, and pointed out to me 
i 4 

] 


that we proceed more systematically by starting with (| 1 
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We have found, then, forms Q'),,, Q?y (the latter real) for which 
all the {a,,,,| are equal (to N-!), B= C= D=1, and |Q| <A in 8. 
The first form exists for every NV, the second for N of the form 2”. 
8. It remains to prove part (2) of Theorem 1. We take a sequence 
of numbers JN, increasing very rapidly with s, and of the form appro- 
priate to QW. or Q2’y, say the former. Writing for brevity 


Qs = Ws, 


oo 
Q* a, Pd as 2, Aan m Yn- 
s= 


m, 


let 


Evidently Q* is bounded in S by an A (since each Q, is). Let us 
denote the range N,_, <m,n<N, of m, n by R,. Then we have 
for m, n of R,, if the increase of N, is rapid enough, 


= lo-Sanin (6-41) Aap. | Bo *lOn|— (0-1) Aa... 


mn mn 
= s-*N-1—(¢+1)-2Nz!,—... > As-2N 
Hence > |a*®,,|* > Asi. 


st 
Rs 
If the increase of N, is rapid enough, we have 
> si (m,n > N,_,) 


Tmn - 
Wists ee A 


> lar 


1 
mni 
Rs 


S la*_ |! a a's 
a [Ann | Tmn Z > i 0; 


m,n s Rs 


and the series on the left diverges. Similar reasoning shows that 
Cm Um Aiverge if the increase of NV, is rapid enough. 


» ) S 
7 Pns 


The systems (e,,,,) of direction-cosines of equal moduli that we 
=e e been led to consider raise a number of interesting problems. Let 
us denote by s,, the maximum value of | > @,,,,| for varying ». How 


~ 


m=1 
small can we make the set of numbers §;, 89, ..., y by choice of 
the system? In the system 
2mnzi 
enn = = N- 
s, is of the order Nin-}, and N!-€ of the s, are large with N. In 


Mr. Ursell’s system (derived from 


e- 
(1) 
by the process explained in § 8), one s,, (viz. s,) is N, two are 3N, 
four are }N, and so on; N!~- are large with JN, as Liehone. Does this, 
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or something like it, always happen? That such questions are not 
without interest appears if we consider the system 
€mn — N~-*x,,(m) 
of § 8 and ignore the principal characters. It is known that 


8, <AlogN; 


n 
more than this is not known, but more can be deduced from a 
generalization of Riemann’s hypothesis concerning the {-function. 


In the opposite direction nothing appears to be known that is not 


trivial ; we do not know that a single s,, is large. 
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